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With  the  trend  toward  higher  shaft  speeds  to  increase  airflow  and  power  output  without 
increasing  physical  size,  there  have  developed  rotor  dynamic  requirements  peculiar  to  the 
class  of  small  turboshaft  engines  being  developed  for  the  Army;  e.g.,  the  short  blades 
require  extremely  small  tip-to  housing  clearances  to  maintain  high  aerodynamic  efficiency. 

This  in  turn  requires  very  small  rotor  shaft  excursions,  to  avoid  blade-housing  interferences, 
which  is  incompatible  with  low  dynamic  bearing  loads  at  high  speeds.  As  shaft  speeds 
increase  it  becomes  more  difficult  to  keep  critical  speeds  (synchronous  whirl)  out  of  the 
operating  range.  In  order  to  operate  safely  through  or  near  the  flexural  critical  speeds,  a * 

rotor  shaft  must  either  be  well  balanced,  highly  damped,  or  both.  Nonsynchronous  whirl 
induced  by  internal  friction  is  probably  the  most  common  type  cf  self-excited  whirl  in 
turboshaft  engines.  The  destabilizing  forces  generated  by  shaft  splines  appear  to  be  one  * 

of  the  most  common  causes  of  this  type  of  response. 

This  report  provides  the  details  of  a program  designed  to  improve  the  analytical  tech- 
niques currently  available  to  Government  and  industry  for  the  design  and  analysis  of 
rotor  systems  incorporating  squeeze  film  bearing  dampers  and  flexible  bearing  supports  to 
control  synchronous  and  nonsynchronous  whirl.  One  product  is  a squeeze  film  bearing 
damper  analysis,  “Damper  Analysis  by  Numerical  Differentiation  and  Integration"  (DANDI), 
capable  of  predicting  pressures  and  response  of  all  known  damper  configurations  under 
normal  and  abusive  unbalance  conditions  and  under  nonsynchronous  response  conditions; 
another  is  a spiine  coupling  analysis  program  that  predicts  the  destabilizing  forces  of 
commonly  used  spline  configurations  with  and  without  lubrication;  and  a third  is  valida- 
tion resting  results  showing  that  existing  beam  formulas  for  calculating  the  stiffness  of 
flexible  bearing  supports  have  satisfactory  accuracy  for  estimating  support  performance  as 
part  of  an  overall  rotor  system  dynamic  analysis. 

Mr.  Allen  C.  Royal,  Propulsion  Technical  Aree,  Aeronautical  Technology  Division,  served 
as  Project  Engineer  for  this  project. 


DISCLAIMERS 

The  findings  in  this  r sport  are  not  to  be  construed  » »n  official  Department  of  the  Army  position  unless  so 
designated  by  other  autf  nrizad  documents. 

When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than  in  connection 
with  a definitely  related  Government  procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  end  the  fact  that  the  Government  may  have  formulated,  furnished, 
or  in  eny  way  supplied  the  said  drawings,  specifications,  or  other  data  is  not  to  be  regarded  by  implication  or 
otherwise  as  in  any  manner  licensing  the  holder  or  eny  other  person  or  corporation,  or  conveying  any  rights  or 
permission,  to  manufacture,  use,  or  sell  eny  pstented  invention  that  may  in  eny  wpy  be  related  thereto. 

Trade  names  cited  in  this  report  do  not  constitute  an  official  endorsement  or  tpproval  of  the  use  of  such 
commercial  hardware  or  software. 


DISPOSITION  INSTRUCTIONS 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to  the  originator. 


wnT7W7Ew*noTXRiw.\*i  "wnwnw.mwin  'W- 


UNCLASSIFIED 


1IFICATIOH  OF  THIS  PAQE  (Wi«n  Dmtm  Snr«« 


fir MOHlto^INQ  AGENCY  NAME  A AODRCSm'/r  dl «•  *•»!  I Inw  Cm/ro ll/nd  Olllct) 


REPORT  DOCUMENTATION  PAGE 


RFAP  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I2.  OOVT  ACCESSION  NO.j  3.  RECIPIENT'S  CATALOG  NUMBER 


fNGINE  ROTOR  DYNAMICS.  I 

YNCHRONOUS  AND  NGNSYNCHRONOUS  \ 

WHIRL  gONTROL.  **  ) 

tT"auThoS7«J 

"I 


armol  J 


.1  hum  «r  nnnrrr~T*y'  ijmj  n n ari 
Final ^Technical  ^fepact  I 


, REPORT  NUMBER 


R GRAN  .'  NUMSCRf*) 

DAAJj0(2-76-C^011  / 


t.  PERFORMING  ORGANIZATION  NAME  ARC  ADURESS 
Pratt  a WhitnFV  Aircraft  Group  / 

Govunmant  Products  Division  v 
P.0.  Bos  2691 

Wsst  Palm  Baach,  Florida  33402  


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a ?ORK  UNIT  LUMBERS 


62C109A  1F2622D9 


AH73  00  144  1^( 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Applied  Technology  Laboratory 
U.  S.  Army  Research  and  Technology 
Laboratories  (AVRADCOM)  Fort.  Eustia,  VA  23604 


A 


MKuariMJJA — , 

Febi-wraynpr/J  I 


S.  NUMBER  OF  PAGES 


148 


'£)  Jr 3 ij 


IS.  SECURITY  CLASS,  fo/ 1W«  raportj 

UNCLASSIFIED 

“OECLAS3IFICATION/ DOWN  GRADING- 


ISa. 


SCHEDULE 


Rotor  Dynamics 
Squeeze  Film 
Damper 
Critical  Speed 


Blade  Loss  Response 
Spline  Coupling 
Stability 

Subcritical  Shaft 


Supercritical  Shaft 
Flexible  Supports 
Nonaynchronous  Whirl 
Synchronous  Whirl 


ABSTRACT  fOrarfM  ■—  r. 


■Mi  tt  mmd  Mmmtttr  Mr  W»ek  niitwj 

, A combined  analytical  and  test  program  has  been  performed  to  develop  a method  of 
designing  high-speed  power  turbine  rotors  to  (1)  minimize  rotor-induced  dynamic  loads  under 
normal  operating  conditions,  (2)  minimize  rotor  tip-to-shroud  clearance  to  maintain  high  flow- 
path  efficiency,  and  (3)  minimize  titor  deflections  due  to  sudden  abusive  imbalance  loads 
associated  with  blade  lots.  A deaign/prediction  system  was  established  for  both  synchronous  and 
nonaynchronous  rotor  whirl  and  the  best  compromise  design  for  a rotor  subjected  to  these 
excitations.  ^ 


■MUM  OF  • NOV  as  IS  OiBOH 


UNCLASSIFIED 


la.  mSTRIRUTION  STATEMENT  (ol  tfll#  RaportJ 

Approved  for  public  release;  distribution  uni 

(fr)  . 

imited. 

17.  DISTRIBUTION  STA1  EMENT  (o  1 thm  mb Jtfmct  mntmtmd  In  Block  20,  IS  dlftmrmttt  from  Rmport)  X'jf 

tt.  SUPPLEMENTARY  NOTES 

SE.CUMTV  CLASSIFICATION  OF  TNI*  PAGE  Dj*f  Snl*r*« 


3/3  33/ 


79  03  19  0" 


T 7JJ& 


UNCLASSIFIED 


...  tv 


:UAl«flCATIO«  or  THIt  P*OI(TWl»>  CM  «nl«f»0 


Two  mathematical  models  were  developed  in  the  course  cf  the  program.  First,  a squeeze 
film  damper  model  and  subsequent  computer  program  were  developed  to  quantify  the  effects  of 
end -seal  leakage  and  inlet  feedback  on  damper  performance.  Secondly,  a spline  coupling  friction 
mathematical  model  and  computer  program  were  developed  to  predict  the  destabilizing  forces 
which  excite  rotor  nonsynchronous  whirl. 


To  verify  the  accuracy  of  the  analytical  models  developed  in  the  program  and  other  existing 
rotor  dynamic  models,  a comprehensive  test  program  was  run.  This  program  included:  (1)  a high 
speed  squeeze  film  damper  test  rig,  (2)  bearing  flexible  support  testing,  (3)  a synchronous 
response  rotor  rig,  and  (4)  a nonsynchronous  response  rotor  rig.  Using  thi  results  of  the  analytical 
models  and  experimental  tests,  a method  of  design  optimization  was  del  doped  to  obtain  the  best 
trade-off  between  all  the  rotor  design  variables  considered  in  this  program. 


IV 


1 


v 


# 


UNCLASSIFIED 


SCCU*lT<r  CL  ICATIO*  or  THIS  Oaf*  Bmt+vm*) 


f!  9*S 


W-jnr.v.- ' 


PREFACE 


This  report  covers  the  development  of  a design/prediction  system  for  optimizing  the  design 
of  high-speed  turboehaft  engines.  The  36-month  Engine  Rotor  Dynamics  Program  ww  performed 
under  Contract  DAAJ02-76-C-0011  with  the  Applied  Technology  Laboratory,  U.  S.  Army 
Research  and  Technology  Laboratories  (AVRADCOM),  Ft.  Eustis,  Virginia.  The  program  was 
performed  under  the  technical  cognizance  of  Mr  A.  Royal  of  the  Propulsior  Croup,  Technology 
Applications  Division,  Applied  Technology  laboratory.  The  results  of  the  program  were 
culminated  in  a method  of  design  optimization  for  power  turbine  rotor  dynamics. 

The  program  was  conducted  by  United  Technologies  Corporation,  Pratt  & Whitney  Aircraft 
Group  Government  Products  Division,  West  Palm  Eeach,  Florida  with  the  subcontracted 
assistance  of  Pratt  & Whitney  Aircraft  of  Canada  LTD.,  Longueuil,  Queoec;  Mechanical 
Technology  Incorporated,  Latham,  New  York;  and  Dr.  J.  M.  Vance,  formerly  Associate  Professor 
of  the  Department  of  Mechanical  Engineering,  University  of  Florida,  Gainesville,  Florida. 
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SECTION  I 


INTRODUCTION 

Recent  trends  in  aircraft  gas  turbine  engine  design  have  been  directed  toward  lightweight, 
flexible  rotor  and  case  systems  in  the  interest  of  increased  horsepower/thrust-to-weight  ratios. 
Without  precise  control  of  rotor  and  case  system  dynamic  deflections,  severe  performance 
deterioration  will  occur  as  a result  of  blade  tip-to-ahroud  rubs.  As  much  as  1%  loss  in  efficiency 
due  to  a surge-induced  rotor  tip  rub  can  occur  in  a multistage  compressor.  More  than  3%  engine 
performance  deterioration  can  result  from  a combination  of  severe  compressor  and  turbine  rotor 
tip-to-shroud  rubs. 

The  rotor  system  dynamic  design  analyses  must  also  address  sudden,  large-scale  rotor 
imbalance.  For  example,  the  loss  of  an  engine  and  possibly  a catastrophic  loss  of  an  aircraft  could 
occur  as  the  result  of  inadvertent  blade  loss,  if  the  rotor/support  system  is  too  sensitive  to 
imbalance. 

The  Applied  Technology  Laboratory  recognized  the  requirement  for  improved  analytical 
design  and  prediction  techniques  in  the  area  of  engine  rotor  dynamics  and  contracted  with  Pratt 
& Whitney  Aircraft  to  conduct  this  comprehensive  development  program.  The  improvement  of 
engine  rotor  dynamics  design  tools  in  this  program  has  three  major  objectives:  (1)  to  minimize 
rotor-induced  dynamic  loads  under  normal  operating  conditions;  (2)  to  minimize  achievable  rotor 
tip-to-shroud  clearance  to  maintain  high  flowpath  efficiency;  and  (3)  to  minimize  rotor 
deflections  due  to  sudden  abnormal  imbalance  loads  associated  with  blade  loss  and/or  engine 
surge. 

The  36-month  program  is  organized  into  10  tasks  (Figure  1).  The  first  seven  tasks  are 
concerned  with  establishing  a design/prediction  system  for  synchronous  whirl  situations.  The  last 
three  teaks  are  directed  toward  nonsynchronous  whir!  and  the  best  compromise  design  for  a rotor 
system  subjected  to  both  types  of  excitation. 

Tasks  I,  II,  and  III  established  the  operating  characteristics  of  an  oil  film  damper.  Task  IV 
established  the  design  characteristics  of  flexible  bearing  supports.  The  combination  of  oil  film 
and  support  structure  characteristics  provided  the  basis  for  a iesign/prediction  system  for  oil- 
damped  rotor  systems.  Task  I,  Analytical  Prediction,  consisted  of  the  formulation  of  an  analytical 
model  of  a squeeze  film  damper.  This  mode!  included  the  Tondl  and  the  Lund  coefficients.  A 
viscous  damper  rig,  which  has  the  capability  to  provide  all  the  experimental  data  needed  to 
determine  the  Tondl  and  Lund  coefficients,  was  designed  in  Task  II,  Conceptual  and  Detail 
Design.  In  Task  III,  Fabrication  and  Test,  bearing  damper  pressures  and  shaft  motion  were 
measured  over  a range  of  design  parameters  that  permitted  data  comparison  with  current  bearing 
theories.  These  data  fully  characterized  viscous  damper  behavior  under  synchronous  loading 
conditions.  In  Task  IV,  establishment  of  the  operating  characteristics  (static  and  dynamic 
stiffness  and  structural  damping)  of  three  different  types  of  bearing  supports  was  accomplished. 
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Th*  validi'v  of  this  equof  se  film  damper  design  system  was  demonstrated  in  an  existing 
con<;rol!od  test  rig  in  Tasks  V,  VI,  and  VII,  which  involved  the  prediction  and  subsequent 
measurement  of  bearing  forces  and  rotor  tc-shroud  deflection.  In  Task  V,  Analytical  Prediction, 
the  Task  I design  system  was  u**d  to  define  an  optimum  damped  bearing  configuration. 
Predictions  were  made  for  dynamic  nearing  forces  and  rotor-to-shroud  displacement  as  a function 
of  six  variables:  (1)  rotor  unbalance,  (2)  rotor  shaft  speed,  (3)  rotor-to-shroud  clearance,  (4) 
bumper  clearance,  (5)  bearing  damper  support  flexibility,  and  (6)  the  most  important  oil  film 
variables  established  during  the  viscous  damper  characterization.  These  predictions  were  made 
for  both  subcritical  and  supercritical  shaft  speeds.  In  Task  VI,  Measurement  Capability,  an 
experimental  rig  was  designed  to  provide  the  capability  to  simultaneously  measure  dynamic 
bearing  loads  and  rctor/shroud  relative  deflections.  The  experimental  rig  was  fabricated  and 
tested  in  Task  VII,  Fabrication  and  Test,  using  the  capability  designed  in  Task  VI  to  measure 
rotor  dynamic  response  as  a function  of  the  variables  specified  in  Task  V.  The  test  results  show 
that  rotor/shroud  deflection  can  be  minimized  by  either  a small  clearance  or  large  clearance 
damper.  However,  the  small  clearance  damper  produces  high  bearing  loads;  therefore  the  large 
clearance  damper  and  bumper  combination  showed  the  best  overall  dynamic  performance.  The 
test  program  also  demonstrated  the  survivability  of  a supercritical  rotor  by  simulating  a blade 
loss  at  a speed  20%  above  the  shaft  bending  critical  speed  and  decelerating  back  again  through 
the  critical  speed  with  the  abusive  blade  loss  imbalance  without  damage  to  the  rig. 

In  Tasks  Vm  and  IX  the  effectiveness  of  the  optimum  damper  system  in  suppressing 
nonsynchronous  excitations  such  as  those  typically  encountered  in  power  turbine  rotor  systems  of 
turboehaft  engines  was  determined.  Task  VIII,  Nonsynchronous  Whirl,  dealt  with  the  problem  of 
predicting  spline-coupling-induced  nonsynchronous  rotor  whirl  and  the  dc  ?ign  of  a test  rig  which 
had  the  capability  to  experimentally  demonstrate  and  measure  nonsynchronous  whirl  and  the 
threshold  speed  of  instability.  A rotor  rig  that  would  demonstrate  nonsynchronous  whirl  as  a 
function  of  the  following  four  variables  was  fabricated  and  tested  under  Task  IX,  Fabrication  and 
Test:  (1)  spline  type,  (2)  spline  clearance  fit,  (3)  spline  lubricant  supply,  and  (4)  damper 
lubricant  supply.  The  test  results  showed  that  squeeze  film  damping  and  snline  lubrication  act 
to  suppress  spline-coupling -induced  rotor  instability.  On  all  occasions  the  test  rig  could  be 
stabilized  by  simply  lubricating  the  spline  teeth  contact  surfaces. 

In  the  final  task,  Task  X,  a trade  study  was  conducted  to  determine  the  best  compromise 
damper  support  design  for  power  turbine  shafts,  considering  both  synchronous  and  non- 
synchronous whirl.  Using  the  results  of  tasks  VI  through  IX,  a method  of  design  optimization  for 
the  beat  trade-off  between  all  the  rotor  design  variables  considered  in  this  program  was 
developed.  Because  of  the  vast  number  of  mission  requirements  for  gas  turbine  engines,  each 
engine  must  be  evaluated  by  the  criteria  dictated  by  its  individual  mission.  However,  a fast 
iterative  procedure  has  been  developed  which  will  allow  rotor  dynamic  optimization  in  the  early 
stages  of  engine  development. 
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SECTION  II 


VI8COU8  DAMPER  CHARACTERIZATION  (TASKS  I,  II,  AND  III) 


A.  GENERAL  DISCUSSION 

A mathematical  model  for  squeeze  Him  dampers  was  developed,  and  the  solution  results 
were  compared  with  data  from  four  different  test  rigs.  A special  feature  of  the  analysis  was  the 
treatment  of  several  different  types  of  end  seals  and  inlets,  with  inlet  feedback  included.  A finite 
difference  method  was  used  to  solve  the  Reynolds  equation,  with  a banded  matrix  inversion 
routine.  The  test  data  were  taken  from  a high-speed  free-rotor  rig  and  from  three  previously 
tested  controlled-orbit  rigs. 

Squeeze  film  dampers  have  been  used  for  many  years  to  minimize  rotor  vibrations  in  gas 
turbine  engines.  The  success  rate  of  these  cLr1  oers  has  been  sporadic,  primarily  due  to  the  lack 
of  adequate  analytical  modeling.  Until  recently,  damper  analyses  were  based  on  a journal  bearing 
theory  that  restricted  the  damper  to  centralized  circular  whirl  orbits  and  simplified  end  flow 
conditions.  The  damper  model  presented  in  this  report  treats  noncircular  orbits,  end  seal  effects, 
and  inlet  feedback  effects. 

In  1886  Reynolds  derived  his  famous  lubrication  equation  based  on  the  theory  that  the 
resultant  pressure  profile  was  due  to  hydrodynamic  action  in  the  fluid  and  was  dependent  upon 
the  viscosity  of  the  lubricant  being  used.  The  Reynolds  equation  was  derived  from  Navier-Stokes 
general  flow  equations  with  assumptions  such  as  thin  film,  laminar  flow,  no  surface  slippage,  etc. 
Several  authors  such  as  Hori  (Reference  1)  and  Tondl  (Reference  2)  analytically  solved  the 
Reynolds  equation  by  assuming  a “long  bearing”  (no  lubricant  end  flow)  for  studying  journal 
bearing  performance.  The  film  circumferential  end  conditions  imposed  by  these  authors  apply  to 
journal  bearings  and  also  to  dampers  with  centralized  circular  whirl  orbits. 

The  computer  program  described  in  this  report  considers  general  damper  orbits  including 
nonzero  radial  velocities  and  the  effects  of  lubricant  flow  out  of  the  end  seals  and  through  the 
inlet.  Piston  ring  and  O-ring  type  end  seals  are  modeled  to  create  continuity  in  flow  throughout 
the  damper,  and  the  resulting  equations  are  written  in  finite  difference  form.  Multiple  hole  and 
circumferential  groove  inlets  are  modeled,  and  the  total  pressure  distribution  is  solved 
simultaneously  by  use  of  a banded  matrix  inversion  technique. 

The  analytical  model  was  verified  by  examining  the  pressure  data  measured  in  four  squeeze 
film  damper  rigs  with  speeds  ranging  to  60,000  rpm.  Three  of  the  rigs  were  previously  reported  in 
References  3,  4,  and  5.  These  three  rigs  were  run  at  low  speeds,  to  3,000  rpm,  and  had  controlled 
orbits  that  precisely  defined  the  motion  of  the  damper.  The  high  speed  data  to  be  described 
herein  were  measured  on  the  fourth  rig  in  which  the  damper  had  unconstrained  orbits  (free  rotor). 
The  four  rigs  have  a wide  variety  of  damper  sizes  and  inlet  and  seal  configurations. 

1 Hori,  Y.,  “A  Theory  of  Oil  Whip,"  Journal  of  Applied  Mechanic*,  Trans.  ASME,  June  1969. 

2 Tondl,  A., “Some  Problems  in  Rotor  Dynamics, "Chapman  and  Hall  Limited,  London,  England,  1966,  pp.  122-130. 

3 Feder,  E„  Banial  P.  N.  and  Blanco,  A.,  “Investigation  of  Squeeze  Film  Damper  Force*  Produced  by  Circular  Centered 
Orbita,"  presented  at  the  Gas  Turbine  Conference,  Philadelphia,  Pa.,  March  27-31, 1977,  ASME  Paper  No.  77-GT-27 

4 Vance,  J.  M.,  and  Kirton,  A.  J.  “Preliminary  Investigation  of  the  Dynamic  Force  Response  Coefficients  for  Squeeze 
Film  Bearing  Dampen,"  U.S.  Army  Research  Office,  Grant  Number  DAHCO4-74-6-0048,  November  1974. 

5 Jonee,  M.  G.,  "An  Experimental  Investigation  of  Squeeze  Film  Hydrodynamics,"  Report  No.  K.  320,  National  Gas 
Turbine  Establishment,  Ministry  of  Defense  (England),  January  1973. 
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B.  ANALYTICAL  MODEL 

1.  Description 

The  problem  ie  that  of  calculating  the  distribution  of  film  pressure  in  a squeeze  film  bearing 
damper,  given  the  journal  position  and  velocity. 

The  approach  ie  to  calculate  the  elemonts  of  the  matrix  and  right  side  vector  that  result 
from  a finite  difference  formulation  of  the  Reynolds  equation  and  the  boundary  conditions.  The 
matrix  and  right  side  vector  define  the  set  of  algebraic  equations  linear  in  the  pressures  at  the 
nodes  of  the  finite  difference  grid. 

A special  feature  of  the  model  is  the  treatment  of  several  realistic  boundary  conditions 
describing  the  inlet  flow  through  any  number  of  holes  or  a groove,  and  the  leakage  flow  through 
several  kinds  of  end  seals  presently  used  in  dampers  for  turbomachinery.  The  types  of  end  seals 
treated  are: 

Piston  t ing 

Radial  O-ring  (around  OD  f journal  at  the  ends) 

Side  nr  axial  O-ring  (again,  the  end  face  of  the  journal) 

No  end  seal. 

Since  the  finite  difference  method  of  solution  requires  no  simplifying  assumptions,  the 
damper  geometry  is  not  restricted;  i.e.,  “short  bearing”  or  “long  bearing”  assumptions  are  not 
required. 

The  solution  efficiency  is  maintained  by  using  a banded  matrix  inversion  routine  (Reference 
S)  to  solve  the  system  of  simultaneous  pressure  equations  resulting  from  the  finite  difference 
breakups.  The  resulting  damper  forces  are  calculated  by  numerically  integrating  the  pressure 
profile. 

2.  Method  of  Solution 

The  squeeze  film  damper  is  considered  to  be  an  oil  film  journal  bearing  with  a nonrotating 
journal.  Thus,  all  of  the  surface  velocity  at  the  film  boundary  is  the  result  of  shaft  translation 
(radial  and  tangential)  allowed  by  the  damper  clearance.  Classical  hydrodynamic  bearing  theory 
shows  that  the  distribution  of  film  pressure  is  given  by  the  solution  of  the  Reynolds  equation, 
subject  to  the  appropriate  boundary  conditions. 


Using  the  notation  and  geometry  from  Figures  2 and  3,  the  Reynolds  equation  is 


6 King,  H.  H.,  "A  Poi*aon  Equation  Solver  for  Rectangular  or  Annular  Regiont,"  International  Journal  for  Numerical  ! 

Method*  in  Engineering,  Vol.  10,  1978,  pp.  799-807.  £ 


Fig-ure  2.  Bearing  Notation  and  Coordinate  System 


Equation  (1)  is  taken  from  equation  f 83)  of  Referenc  : 7,  with  the  velocity  tj.,  taken  as  zero 
(at  the  stationary  bearing  surface).  The  equation  is  made  diraensionlees  by  substitution  of  the 
following  dimensionless  parameters: 


, _ X 7 

v -tc2 
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6/i  R|  iii 
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rj  = _A_,  o = -(^- 
o>R*  ' * wC 

The  derivatives  of  H,  H*.  and  0,  with  respect  to  \p  *re  evaluated  in  functional  form. 
Equation  (1)  then  becomes 
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7 TrumpJer,  P.  R.“Deai(n  of  film  BMrinp!*MacMiiUn  Co.,  New  Yak,  1968,  p.  30. 


and 


T = ut,  t = time  (sec). 


The  damper  pressure  distribution  is  the  solution  to  equation  (2),  subject  to  the  boundary 
conditions  as  determined  by  the  inlet  and  by  the  seals. 

3.  End  Seal  end  Inlet  Boundary  Condition* 

Piston  Ring  Seal.  Figure  4 shows  the  configuration  of  this  seal. 


The  axial  flow  balance  across  any  unit  circumferential  length  around  the  seal  is  given  by 

(3) 


h*  i£_  - r <P-P«)h* 

12/u  aZ  Cp  12mw  ,n/sec/in- 


where  p and  ap/aZ  are  the  pressure  gradient  at  the  inside  face  of  the  seal  (any  X,  Z = L);  p*  is  the 
pressure  outside  the  seal. 

In  dimensionless  form,  this  boundary  condition  is  written 

(*•£)- c’  [p  ( *•£)  -p-]  • 


where 


CP'  = 


. «A» 

p wh* 


(4) 


Note  that  Cp  varies  around  the  seal  with  1/h*. 

Radial  O-ring  seal.  Figure  5 shows  this  seal  configuration. 

Since  the  O-ring  is  flexible,  it  deforms  in  response  to  the  local  film  pressure.  A traveling 
wave  of  axial  deflection  follows  the  rotating  pressure  distribution  and  accommodates  axial  flow 
from  the  film,  even  if  no  leakage  occurs  (Figure  6). 
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(A)  TRAVELLING  WAVE 


.10  ■ h0R 


i(X)  - kor  p(.X) 

*■  <x>  “ K0R  j£ 

(6)  VOLUME  DISPLACED  AT  POINT  A IN  TIME  .IT 

Figure  6.  Traveling  Wave  Around  Radial  0-Ring 


Including  leakage,  the  axial  flow  balance  across  any  unit  circumferential  length  around  the 
sea!  is 


- ho«  K0„  0 R*  -g-  t'X,  L)  - 1— - (X,  L) 


= C0R  [p  (X,  L)  - pe]  . (5) 


In  dimensionless  form,  this  boundary  condition  is  written 


where 


and 


r>  - _ r 

— V-'OR  > 


, 12/i  R*  H0r  0 T, 

~ jjj  e.-oR 


Note  that  CR'  and  K„'  vary  around  the  seal  with  1/hf  and  that  K„'  varies  with  the 
instantaneous  orbit  velocity  0. 

Side  O-ring  seal.  Figure  7 showB  this  type  of  seal.  Notice  that  the  dimensions  h,  and  “a” 
form  a circumferential  cavity  at  the  end  (Z  = L)  of  the  film. 
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This  boundary  condition  requires  that  the  flow  cut  of  the  cavity  (across  the  O-ring)  must 
equal  the  integral  around  the  circumference  of  the  axial  flow  from  the  film  at  Z = L.  Since  h,  » 
hj,  there  is  no  circumferential  pressure  gradient  around  the  cavity.  The  flow  balance  is 

- -jg-  (X,  L)  Rt,  dX  = Cc  [pc  - p.l,  (7) 

with 

p (X,  L)  = p.  , 

where  p and  dp/sZ  are  the  pressure  and  axial  pressure  gradient  at 
Z = L,  ( = 0 at  Z = L ) . 


In  dimensionless  form,  this  boundary  condition  becomes 
- f (#.$  <V-C.  [P.  P.l 


with 


U-fe-)- 


Pc.  all 


where 


C ' = C 

c* 


(8) 


No  end  seal.  If  the  damper  end  is  open,  the  boundary  condition  is 

P (X,  L)  = p.  . (9) 


In  dimensionless  form,  this  becomes 

p ( *■  it ) - p-  • 


(101 


Single  hole  with  backflow. 
The  flow  balance  is 
Qi  ■ Q.  inf/sec, 
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where 


Q,  = flow  into  film  from  inlet, 

Q,  = flow  away  from  inlet  location. 


In  terms  of  pressures  and  gradients,  thio  boundary  condition  is 
Qi  = C,  [p,  — p (X|,  Z|)  ] 

-0--T5T  [-S'  (x, -4-’z'  ) 
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r ap 
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I ] dZ 

(11) 

In  dimensionless  form,  with  d,  < A,  this  becomes 

-f-  w,  ?,)  i.  - -f - »•  +.  i.)  a. 

*»p  up 

+ (4>i,  2,  -)  A^,  — — (^i,  Z,  +)  A^ 

= cr  [p.  - p(*„  ?,)  i 


where 


C ' = - r n 

jjt  '■o  > 


Equation  (12)  applies  at  each  inlet  hole  location. 
Circumferential  inlet  groove  (360°).  The  boundary  condition  is 
p (X,  Zi)  = p,  . 

In  dimensionless  form,  this  becomes 

P (*.  2.)  = P.  • 


PPPUf^g  P^PPllpWSBq8IWWPiPP|r|)^^ 


C.  COMPUTER  PROGRAM 

1.  Flntt*  Difference  Form 

The  computer  program  calculates  the  elements  of  the  matrix  and  right  side  vector  from 
three  sources: 

1.  The  dimensionless  Reynolds  equation  (2),  converted  to  finite  difference  form. 

2.  One  of  the  dimensionless  boundary  conditions  (4,  6,  8,  or  10)  for  the  chosen 
type  of  end  seal  (for  each  of  the  two  damper  ends),  converted  to  finite 
difference  form. 

3.  One  of  the  dimensionless  boundary  conditions  (12  or  J4)  for  the  chosen  type 
of  inlet,  converted  to  finite  difference  form. 


The  equations  indicated  above  are  converted  to  finite  difference  form  by  substitution  of  the 
following  difference  expressions  for  the  derivatives: 


(15) 

(16) 

(17) 


(18) 


where  P;k  is  the  dimensionless  pressure  at  the  node  j,k  of  the  finite  difference  grid,  = 
(j-l)A^,  Z=  (k-l)A„  and  A is  the  grid  size. 

The  grid  size  is  determined  by  choosing  the  number  of  points  M,  N to  be  defined  on  the 
damper  circumference  and  length,  respectively. 

= JFT  ’ = TTT  "rT  (l9) 

A concise  description  of  the  finite-difference  method  as  applied  to  hydrodynamic  bearing 
theory  can  be  found  in  Reference  8. 


8 Pirkui,  O.,  and  Stemlicht,  B.VTbeory  of  Hydrodynamic  Lubrication! 'McGraw-Hill,  Inc.,  New  York,  1961,  pp.  80-81. 
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The  resulting  finite-difference  equations  are  as  follows: 

Main  Program 

Pj*k-X  + [ 1 ~ flj.k  ]PiU  - 2 [ -|f-  + 1 ] PjJ 


+ 


"if  Pj-k+1  + [ 1 + fli.k  ] Pi+1'k  - A/Rj., 


(20) 


where 


^,k  f3jk  dT  + f4i.k  dT  ’ 


and 


f3j.k  - ^j.k 

are  the  functions  associated  with  equation  (2),  evaluated  at  the  point  j,k. 
Subroutine  Piston  (Piston  ring  seai: 

Pj.N  [1  + A.Cp  l -PjiN.i  = AtCpP. , 

where 


N = 


L 

RbA, 


+ 1. 


Subroutine  Rador  (Radial  0 -ring  seal) 


KBpj.i,N 


(22) 


where 
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Subroutine  Sidor  (Side  O-ring  seal) 

M ^ 

L {H'j.N-l  (Pc  - Pj,N-l  )}  -=*-  + CcPc  = Cr  Pe  , (23) 

Az 

J-l 


with 


Pj,N  = Pc  . 

where 


N = 


L 

R*A, 


+ 1. 


Subroutine  Noseal  (No  end  seal) 


fj.N 


= P. 


(24) 


(25) 


where 


N =• 


R*A. 


■ + 1. 


Subroutine  Holes  (Inlet  holes  with  backflow) 


Pj„Nl  [ 2 + 2 + "jf"  <V  J - Pjrl.NI  ~ Pjp  + 1,N1 


A^ 

Pj„NM  “ — 


A,* 


A,1 


A* 

Pji,Nl  + l = C,  P.  , 

A, 


where 


C ■ = J2/i  c, 

' h’  ' ’ 


j = J,  defines  the  circumferential  location  o.  the  ith  inlet 
hole. 

k = N i defines  the  axial  location  of  all  the  inlet  holes. 


(26) 


The  elements  of  the  coefficient  matrix  A,  the  vector  RHS,  and  the  solution  vector  P are 
shown  in  equation  (27).  The  banding  of  the  coefficient  matrix  is  accomplished  by  numbering  the 
grid  points  of  the  finite  difference  grid  in  an  alternating  fashion  around  the  circumference  of  the 
damper. 


MATRIX  EQUATION  FOR  M-6,  N«4 
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f A]  | P | = | RHS  | 

(Matrix  coefficients  are  numbered  as  single  array) 

2.  Solution  of  Proaeure  Matrix 

The  solution  algorithm  is  a Gaussian  elimination  of  the  elements  below  the  diagonal  and  a 
back  substitution  of  unknowns  in  the  reverse  order.  Integration  of  the  pressure  profile  to  get 
damper  forces  is  done  by  Simpson’s  Rule  for  an  even  number  of  finite  strips. 

Premures  that  fall  below  the  damper  cavitation  pressure  are  set  equal  to  the  cavitation 
premure  for  calculating  the  damper  forces.  Film  cavitation  is  not  well  defined;  however,  test 
results  show  that  cavitation  premure  is  approximately  - 14.7  peig  (Reference  3). 

The  numerical  accuracy  of  the  above  procedure  was  verified  by  comparing  the  predicted 
premures  with  known  "long  bearing”  and  "short  bearing”  solutions  to  within  2%.  The  finite 
difference  model  had  a grid  of  6 X 40  end  was  analyzed  in  double  precision  arithmetic  on  an  IBM 
370  Model  168  computer.  The  average  computation  run  time  for  this  size  model  is  approximately 
0.1  sec. 
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A description  of  the  complete  computer  program  (Report  No.  PR- 10732)  waa  provided  to  the 
Applied  Technology  Laboratory  in  October  1978.  A users  manual  for  this  program  appears  in 
Appendix  A and  B. 

D.  EXPERIMENTAL  RIO  DESIGN 

The  objective  of  the  high-speed  squeeze  film  damper  test  was  to  measure  bearing  damper 
forces  and  pressures  over  a range  of  design  variables  to  enable  comparison  of  the  data  with 
existing  theories  for  predicting  damper  performance.  The  test  rig  duplicated  the  dynamics  of 
existing  engine  squeeze  film  dampers,  including  the  case  of  nonzero  radial  velocities. 

1.  Damper  Design 

The  squeeze  film  damper  rig  configuration  is  shown  in  Figure  8.  Design  modifications 
incorporated  into  the  existing  damper  rig  (Reference  9)  are  summarized  in  Table  1.  These 
modifications  permitted  the  measurement  of  damper  squeeze  film  pressures  as  a function  of  shaft 
position  and  velocity,  and  offset  the  shaft  orbit  to  produce  nonzero  radial  velocities. 

The  test  rig  consists  of  two  accurately  balanced  disks  joined  together  by  a hollow  shaft 
supported  at  the  middle  by  a double  row  ball  bearing.  This  mounts  to  the  inner  damper  ring.  The 
disk  assembly  is  driven  by  a 2-in.  air  turbine  via  a quill  shaft  with  a spline  fitting  on  the  driving 
end.  To  induce  motion  into  the  damper,  the  disks  are  unbalanced  by  a known  amount.  The 
damper  force,  pressure,  and  displacement  measurements  are  then  taken  at  various  rotational 
speeds  while  the  disks  spin  in  an  evacuated  chamber. 

The  outer  race  of  the  duplex  ball  bearing  is  shrunk-fit  inside  the  damper  journal  which  is 
held  in  the  damper  bearing  by  a retaining  ring  and  hac  complete  lateral  freedom  of  movement. 
The  damper  journal  has  a nominal  outside  diameter  of  7.6  cm  (2.99  in.).  The  damper  ^as  a length 
of  2.5  cm  (0.98  in.)  divided  equally  between  the  top  and  the  bottom  bands  with  a central  oil- 
groove  of  about  1.25  cm  (0.49  in.)  width.  This  rather  wide  groove  was  required  to  fit  the  off- 
centered  oil  supply  holes  while  retaining  a symmetric  damper.  The  oil  film  is  limited  in  length  at 
each  damper  end  by  a steel  piston  seal  ring. 

A radial  spring  arrangement  is  provided  to  simulate  the  gravity  effect  that  exists  in  a 
horizontal  rotor.  The  side  load  on  the  journal  can  be  controlled  by  adjusting  the  length  of  the 
radial  springs  to  produce  a resultant  pull  of  5.6  kg  (12.4)  on  the  damper,  equal  to  the  wei  of 
the  rotating  parts. 

The  clearance  space  between  the  damper  sleeve  and  the  bearing  housing  is  filled  with  the 
damping  fluid  which  is  pressure  fed  into  the  central  circumferential  groove. 

The  rotor  is  mounted  vertically  in  a containment  casting  that  is  evacuated  to  60  Torricil. 
The  rotor  is  driven  by  a commercial  5-cm  (2-in.)  air  turbine  via  a quill  shaft.  The  rotor  can  be 
driven  over  a wide  speed  range.  The  maximum  speed  is  limited  by  the  bearing  temperature  and 
did  not  exceed  50,000  rpm  during  the  tests.  The  excitation  force  is  provided  by  two  equal  weights 
attached  to  the  top  and  bottom  disks  in  the  same  vertical  plane. 


9 Botraan,  M.,  "Experiment!  on  Cil-Film  Dampen  for  Turbomachinery,”  Journal  of  Engineering  for  Power,  Trans. 
ASME.  Series  A.  Vot.  96.  July  1976,  pp.  393-399. 
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Figure  8.  Squeeze  Film  Damper  Rig 


2.  Inetru  mentation 

The  oil  supply  line  to  the  damper  is  equipped  with  a regulator  valve,  a pressure  gage,  and 
a flowmeter.  The  inlet  oil  pressure  can  be  adjusted  in  the  range  of  5 to  100  peig.  The  oil  coming 
out  of  the  damper  paaeee  through  another  flowmeter.  The  two  flowmeters  provide  information  on 
the  flowrate  within  the  damper  as  well  as  the  leakage  from  its  seals. 

The  oil  film  and  the  bearing  temperatures  are  monitored  by  two  ccppcr-con  thermocouples 
mounted  flush  with  the  damper  journal. 

The  transmitted  load  is  measured  by  the  two  strain  gage  bridges  placed  at  right  angles  to 
each  other  on  the  curved  spokes  of  the  load  cell.  The  load  cell  was  calibrated  statically  to  measure 
radial  loads  in  two  directions.  It  has  a radial  stiffhess  of  243  x 10*  N/m  (1.39  X 10*  lb/in.).  The 
output  signal  from  the  strain  gages  is  amplified  and  record sd  on  an  FM  magnetic  tape. 
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TABLE  1 

DAMPER  RIG  MODIFICATIONS 


Item 

Revision 

Damper  Sleeve 

Provision  was  mode  for  offset 
orbits  and  circumferential  oil 
groove. 

Pressure  Probe  Bushing 

Will  now  accept  dynamic 
(Kistler)  and  static  (Statham) 
pressure  probes. 

Load  Cell 

Additional  pressure  probe  taps 
were  included. 

Load  Cell  Support 

Additional  holes  for  securing 
spring  device  which  produces  off- 
set orbits. 

Retainer  Ring 

Ring  material  was  changed  from 
steel  to  Teflon  to  reduce  friction. 

Dumper  Seals 

Seal  material  was  changed  from 
Teflon  to  steel  to  be  more  repre- 
sentative of  typical  damper  de- 
signs. 

The  displacement  of  the  damper  sleeve  at  the  oil  film  is  measured  by  two  inductance  type 
proximity  probes  mounted  in  the  load  cell  at  the  midplane  of  the  top  land,  displaced  90  deg  from 
each  other.  Both  ac  and  dc  components  of  these  probes  can  be  monitored. 

Eight  holes  with  a diameter  of  0.126  cm  (0.060  in.)  were  drilled  46  deg  apart  in  the  midplane 
of  the  bottom  land  of  the  damper.  These  holes  are  connected  to  the  Kistler  and  the  Statham  type 
pressure  transducers.  The  Kistler  transducers  are  used  to  measure  the  dynamic  pressure,  and  the 
Statham  transducers  the  static  oil-film  pressures.  Thus,  the  instantaneous  fluid-film  pressure 
distribution,  which  is  rotating  with  the  speed  of  the  rotor,  is  measured  in  the  midplane.  Because 
of  space  limitations,  only  four  holes  were  drilled  in  the  top  land  of  the  damper  for  pressure 
measurements.  The  output  signals  from  the  pressure  transducers  are  recorded  on  magnetic  tape. 

The  inertia  effects  of  the  support  mass  are  measured  by  means  of  two  B&K  accelerometers 
mounted  orthogonally  at  the  load  cell.  The  accelerometer  output  is  also  recorded  on  the  magnetic 
tape  via  charge  amplifiers. 

The  amplitudes  of  the  lower  and  the  upper  disks  are  measured  by  two  inductance  type 
displacement  probes.  The  reference  signal  from  a marker  on  the  rotor  is  used  to  measure  the 
phase  angles  between  different  signals.  A magnetic  probe  is  used  for  this  purpose.  Another 
magnetic  probe  is  used  for  monitoring  the  speed  of  the  rotor. 

The  damper  proximity  probes  were  calibrated  before  each  build  to  study  the  drift  in  the 
calibration  constants  from  build  to  build.  No  significant  drift  was  observed. 
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The  damper  oil  usad  is  a synthetic  oil  MIL-L-23699  whose  viscosity  and  specific  gravity 
varies  as  a function  of  its  temperature.  In  the  tests,  the  temperatures  varied  from  27*to  93°C  (80* 
to  200°F  ; with  the  corresponding  variation  in  absolute  viscosity  from  37.85  x lO-’  N-sec/m*  (5.47 
x 10'f  lb  sec/in.*)  to  6.06  x 10"*  N-sec/m*  (0.73  X 10~*  lb  sec/in.*). 

The  basic  damper  instrumentation  is  shown  in  Figure  9. 


Figure  9.  Damper  Rig  Instrumentation 

E.  RIQ  TEST  AND  RESULTS 
1.  Test  Program 

Tests  were  performed  for  three  values  of  radial  clearance  and  three  values  of  imbalance  for 
both  concentric  and  eccentric  orbits.  Table  2 summarizes  the  parameter  combinations  for  the 
various  runs.  The  effect  of  inlet  oil  pressure  was  also  studied  at  three  different  values  of  the 
pressure. 

The  rotor  was  run  over  a wide  speed  range  and  measurements  were  taken  at  selected 
constant  speeds  at  intervals  of  5000  rpm  up  to  a maximum  speed  of  50,000  rpm.  At  each  speed, 
measurements  were  taken  of  the  dampe.  sleeve  deflections,  transmitted  loads,  support 
accelerations,  instantaneous  oil-film  pressures  and  the  associated  phase  angles  The  oil  flow  and 
bearing  temperatures  were  recorded  manually  at  each  of  these  speeds. 

The  damper  sleeves  had  a circumferential  oil  inlet  groove  and  piston  ring  end  seals  which 
simulated  a typical  high  speed  engine  damper  configuration.  The  damper  radial  clearances  were 
0.006,  0.010,  and  0.015  in.,  typical  of  current  damper  designs. 


TABLE  2 

THE  PARAMETER  COMBINATIONS  FOR  VARIOUS  RUNS 


Run  No. 

Orbit 

Inlet  Oil 
Pressure 

. (psig) 

Imbalance 
(gm  cm) 

Radial 

Clearance 

(mm) 

Imbalance  Phase 
Relative  to  Marker 
(deg) 

1.1 

Concentric 

30 

0 

0.127 

15  lag 

1.2 

Concentric 

30 

27.4 

0.127 

15  lag 

1.3 

Concentric 

30 

13.7 

0.127 

15  lag 

1.4 

Eccentric 

30 

13.7 

0.127 

15  lag 

1.5 

Eccentric 

30 

27.4 

0.127 

15  lag 

1.6 

Eccentric 

30 

0 

0.127 

15  lag 

1.7 

Concentric 

30 

0 

0.127 

15  lag 

2.1 

Concentric 

30 

0 

0.254 

15  lead 

2.2 

Concentric 

30 

13.7 

0.254 

15  lead 

2.3 

Concentric 

30 

27.4 

0.254 

15  lead 

2.4 

Eccentric 

30 

27.4 

0.254 

15  lead 

2.5 

Eccentric- 

30 

0 

0.254 

15  lead 

2.6 

Eccentric 

30 

13.7 

0.254 

15  lead 

2.7 

Eccentric 

60 

13.7 

0.254 

15  lead 

2.8 

Eccentric 

100 

13.7 

0.254 

15  lead 

3.1 

Concentric 

30 

0 

0.381 

15  lead 

3.2 

Concentric 

30 

13.7 

0.381 

15  lead 

2.3 

Concentric 

30 

27.4 

0.381 

15  lead 

3.4 

Eccentric 

30 

27.4 

0.381 

15  lead 

3.5 

Eccentric 

30 

0 

0.381 

15  lead 

3.6 

Eccentric 

30 

13.7 

0.381 

15  lead 

3.7 

Eccentric 

5 

13.7 

0.381 

15  lead 

2.  Test  Results 

Some  of  the  results  obtained  from  the  tests  are  prjsented  in  Figures  10  through  16. 
Figure  10  shows  the  variation  of  the  bearing  and  oil-film  temperatures  with  speed.  The  influence 
of  the  unbalance  on  the  bearing  temperatures  was  found  to  be  marginal,  the  temperatures  being 
slightly  higher  for  higher  unbalances.  The  oil  leakage  through  the  piston  rings  is  plotted  in 
Figure  11.  It  is  suspected  that  most  of  the  oil  leakage  occurred  through  the  location  where  the 
piston  ring  ends  butt  together. 

The  damper  orbits  remained  nearly  circular  and  the  pressure  profile  nearly  sinusoidal  up  to 
a certain  speed  range.  As  the  speed  was  further  increased,  the  pressure  profile  exhibited  a flat 
portion  indicating  a region  with  zero  pressure.  This  is  referred  to  as  onset  of  cavitation.  When  the 
speed  was  further  increased,  the  orbits  became  irregular,  indicating  nonsynchronous  whirl. 
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Oil  Supply  Pressure  (psi) 

Figure  16.  Effect  of  Oil  Supply  Pressure  on  Transmissibility  and  Damper  Deflections 


The  cavitation  was  observed  to  progress  rapidly  and  was  followed  by  multi-loop  orbits.  At 
the  onset  of  multi-loop  orbits,  the  film  extent  was  generally  found  to  be  less  than  180  deg.  The 
multi-loop  orbits  show  up  as  a fractional  frequency  component.  Spectrum  analysis  of  the  signals 
indicated  this  fractional  frequency  to  be  centered  around  270  ± 26  Hz.  This  is  shown  in 
Figure  12,  where  frequency  of  the  nonsynchronous  component  is  plotted  as  a percentage  of 
synchronous  frequency.  After  the  onset  of  this  instability,  the  fractional  frequency  remained 
practically  constant  until  it  reached  half  the  rotor  speed.  At  higher  speeds,  this  frequency  then 
remained  at  half  the  rotor  speed  with  substantially  reduced  associated  amplitudes.  The  ratio  of 
nonsynchronous  vibration  amplitude  to  the  synchronous  vibration  amplitude  is  also  shown  in 
Figure  12.  It  decreased  with  increased  speeds. 

With  a balanced  rotor,  no  fractional  frequency  whirl  was  observed.  With  a moderate  amount 
of  imbalance,  the  fractional  frequency  whirl  disappeared  beyond  a certain  speed  range.  For  all 
rune  with  an  unbalance  of  13.7  gm  cm  (0.19  oz  in.),  the  fractional  frequency  component 
disappeared  at  rotor  speeds  exceeding  43,000  rpra.  The  orbits  became  circular  again  when  the 
fractional  frequency  component  disappeared  at  higher  speeds. 

Figures  13  and  14  show  the  size  and  the  location  of  the  damper  orbit  within  the  clearance 
circle  for  runs  1.3  and  1.4,  respectively,  at  different  speeds.  With  increasing  speeds  and 
unbalances,  the  orbit  size  grew  and  the  journal  center  moved  toward  the  bearing  center.  The 
orbits  tended  to  be  concentric  at  higher  speeds  and  higher  unbalances.  The  orbit  size  and  shape 
for  tho  rotor  with  radial  springs  were  similar  to  those  without  radial  springs  except  that  the 
journal  center  was  displaced  toward  the  redial  spring  force  for  eccentric  orbits.  At  low  speeds,  the 
journal  was  displaced  on  one  side  by  the  impact  of  < he  inlet  oil  jet. 

The  transmissibility,  as  defined  by  the  retio  of  the  load  transmitted  to  the  load  cell  to  the 
unbalance  force,  is  shown  in  Figure  16.  It  is  seen  to  be  lower  for  dampers  with  larger  clearances 
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and  slightly  higher  for  eccentric  orbits  up  to  the  speeds  where  the  oil  film  forces  were  greeter  than 
radial  spring  forces. 

The  effoct  of  inlet  oil  pressure  on  the  transmissibility  and  the  damper  deflections  is  shown 
in  Figure  16.  The  high  inlet  oil  pressure  reduced  both  the  damper  deflections  and  transmitted 
loads  in  the  synchronous  operating  range  of  the  damper  as  shown  in  the  figure. 

F.  DATA  CORRELATION  WITH  ANALYTICAL  MODEL 

The  objective  of  the  data  correlation  was  to  verify  the  ability  of  the  computer  analysis  model 
to  predict  the  forces  and  pressures  that  occur  in  damper  designs  commonly  used  in  gas  turbine 
engines. 

The  test  data  were  taken  from  the  high-speed  free-rotor  rig  test  (Section  II.D)  plus  three 
previously  reported  low-speed  controlled-orbit  rig  tests  (References  3.  4,  and  5).  The  free  rotor 
damper  rig  attained  speeds  of  50,000  rpm,  while  the  maximum  speed  of  the  othar  three  rigs 
(Figure  17)  was  3000  rpm.  These  rigs  represented  a goed  cross  section  of  dampers  commonly  used 
in  the  gas  turbine  engine  industry  by  exhibiting  a variety  of  damper  sizes  and  enri-seal 
configurations.  The  resultc  of  the  computer  analysis  described  above  agree-  closely  with  the  data 
from  the  three  controlled- orbit  rigs. 

1.  Hlgli-Sp?cd  Rig 

Figure  18  shows  the  squeeze  film  pressure  correlation  with  the  high-speed  free-rotor  rig.  The 
damper  speed  is  20,000  rpm  and  the  radial  clearance  is  5 mils.  The  data  correlation  shown  is 
typical  of  all  rig  speeds  and  clearance.  The  eight  measured  data  points  were  recorded  at  an 
instant  in  time  by  the  eight  equally  spaced  pressuie  transducers  placed  around  the  circumference 
of  the  damper.  The  corresponding  predicted  pressures  were  calculated  for  the  same  instant  of 
time.  Tfc?  triangle  symbols  represent  the  pressures  predicted  by  modeling  the  piston  ring  seals 
with  a radial  gap  of  1 mil.  These  pressures  were  always  predicted  higher  than  the  measured 
values,  which  suggests  that  this  type  of  eal  may  leak  more  than  had  been  anticipated.  If  the 
piston  ring  end  seals  are  modeled  as  being  totally  ineffective  (unrestricted  end  leakage),  the 
predicted  and  measured  pressure  compare  very  closely.  The  unrestricted  leakage  assumption  is 
supported  by  the  relatively  high  leakage  rate  of  0.3  gal/min  measured  at  this  test  point. 

2.  Low-Speed  Rigs 

Figures  19  through  21  show  the  comparison  of  measured  and  analytically  predicted  squeeze 
film  pressures.  These  figures  also  show  the  predictions  of  two  currently  used  analytical  models, 
based  on  simplified  assumptions. 

Note  that,  although  the  "long  bearing”  theory  is  ;n  reasonable  agreement  with  one  rig  and 
the  “short  bearing”  tneory  with  another,  the  finite  difference  model  is  in  good  agreement  for  all 
three  controlled -orbit  rigs. 
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8.  REF  4 RIG-SINGLE  HOLE  INLET,  AXIAL  O-RING  SEALS 
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C.  REF  5 RIG-CIRCUMFERENTIAL  GROOVE  INLET. 
SMALL  CLEARANCE  SEAL 


Figure  17.  Low-Speed,  Contmlkd-Crbit  Damper  Riga 
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RADIAL  VELOCITY  - -0.931  In./sec 
TANGENTIAL  VELOCITY  - 0.778  In. /sec 
Damper  (Fig  7) 

RADIUS  - 1.5  In.  (3.81  cm) 

LENGTH  - 1.0  In.  (2.54  cm) 

CLEARANCE  - 0.0052  In.  (0.0132  cm) 
ECCENTRICITY  - 0 0031  In.  (0.0079  cm) 
SUPPLY  PRESS-  30  pslg 

SPEED  - 20,000  rpm 

VISCOSITY  - 3.3X10'°  Ib-sec/ln.2  (reyn) 


Figure  18.  Free -Rotor  Rig  Data  Correlation 

The  controlled-orbit  rig  reported  by  Feder  (Reference  3)  had  O-ring  side  seals  that  allowed 
very  little  leakage.  Therefore,  the  "long  bearing”  approximation  to  the  lubrication  theory  shows 
nearly  as  gooo  a comparison  with  the  experimental  data  as  the  finite  difference  analysis  (Figure 
19).  The  rig  used  by  Vance  and  Kirton  (Reference  4)  also  had  very  effective  end  seals  and  had  a 
single  hole  inlet,  so  it  also  is  approximated  fairly  well  by  the  “long  bearing”  model  as  shown  >n 
Figure  20.  However,  the  predicted  peak  pressure  deviates  from  the  measured  amplitude  by  a 
factor  of  two,  a deviation  which  is  significantly  reduced  by  the  predictions  of  the  finite  difference 
model. 

At  the  other  extreme,  the  rig  used  by  Jones  (Reference  5)  had  a circumferential  oil  feed  and 
end  seals  which  allowed  considerable  leakage.  As  shown  in  Figure  21,  the  “long  bearing” 
approximation  is  completely  inappropriate  for  this  configuration.  The  “short  bearing”  model  is 
a closer  approximation  but  does  not  fit  the  data  as  well  as  the  finite-difference  model  that 
includes  both  the  effects  of  the  end  leakage  and  the  inle*  fiow. 
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ANGLE  MEASURED  FROM  MINIMUM  FILM  THICKNESS  - deg 

DAMPER  1 (Fig  8a) 

RADIUS  - 3.22  In.  (8.18  cm) 

LENGTH  =*  2.0  In.  (S.08  cm) 

CLEARANCE  - 0.020  In.  (0.0508  cm) 
ECCENTRICITY  - 0.014  In.  (0.0356  cm) 

SUPPLY  PRESS-  30  pslg 

SPEED  - 1000  rpm 

VISCOSITY  - 2. 62X10"®  Ib-aec/ln.2  (reyn) 


'igure  19. 


Controlled- Orbit  Rig  Data  (Ref  3)  Correlation 


ANGLE  MEASURED  FROM  MAXIMUM  FILM  THICKNESS  - DFG 


CAMPER  2 (Fig  8b) 

RADIUS 

LENGTH 

CLEARANCE 

ECCENTRICITY 

SUPPLY  PRESSURE 

SPEED 

VISCOSITY 


- 2 5 in  (6  35  cm) 

= 4 813  in.  (12,23  cm) 

« 0 062  In.  ;0  157  cm) 

« 0.031  in  (0  0787  cm) 

« 60  pslg 
* 1766  i pm 

- 27  4X10'®  lb-sec/ln.2  (reyn) 


Figure  20.  Controlled -Orbit  Rig  Data  (Ref  4)  Correlation 
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Angle  Measured  From  Maximum  Film  Thickness  - deg 


Damper  (Figure  17,  Part  C) 


Radius 

Length 

Clearance 

Eccentrlcty 

Supply  Press 

Speed 

Viscosity 


- 1,5  In.  (3.81  cm) 

* 2.0  In.  (5.08  cm) 

* 0.010  In.  (0.254  cm) 

- 0.0068  In.  (C.0173  cm) 

3 10  pslg 

» 3000  rpm 

■ 12.3  x 10‘®  lb-sec/ln.2  (reyn) 


Figure  21.  Controlled-Orbit  Rig  Data  (Ref  5j  Correlation 


In  summary,  the  predictions  of  the  computer  analysis  model  are  in  close  agreement  with 
experimental  data  taken  from  the  three  controlled-orbit  rigs,  and  are  in  fair  agreement  with  data 
from  the  high-speed  free-rotor  rig.  Measurements  taken  from  the  latter  type  rig  are  inherently  less 
precise,*  although  they  simulate  actual  gas  turbine  conditions  much  more  closely.  Overall,  the 
data  correlations  demonstrate  a high  order  of  accuracy  and  versatility  for  the  finite  difference 
model  described  previously. 


'There  are  more  variable*  to  be  measured  in  this  type  rift,  under  more  difficult  conditions. 
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SECTION  III 


BEARING  SUPPORT  CHARACTERIZATION  (TASK  IV) 


Impedance  tests,  using  a shaker  apparatus,  were  used  to  determine  the  stiffnesses,  damping 
values,  and  natural  frequencies  of  nine  flexible  bearing  supports  of  three  different  designs:  round 
rod  supports,  squirrel  cage  supports,  and  curved  beam  supports.  Three  variations  of  each  support 
type  ware  tested.  The  round  rod  and  squirrel  cage  designs  were  chosen  because  of  their 
widespread  use  in  turboshaft  engines;  the  curved  beam  supports,  an  innovative  new  design,  were 
chosen  because  of  their  compactness,  low  stress  levels  under  loading,  and  the  high  amount  of 
damping  they  provide.  The  main  purpose  of  the  experiment  was  to  test  the  validity  of  current 
analytical  techniques  used  in  predicting  support  stiffness;  the  predicted  and  experimental 
stiffness  data  correlated  well.  A formula  based  on  beam  theory  was  used  in  predicting  the  stiffness 
of  the  round  rod  and  squirrel  cage  supports,  and  Castigliano’s  Theorem  was  combined  with 
Curved  Beam  Theory  to  yield  the  stiffness  of  the  curved  beam  supports.  The  experimental  results 
and  the  analytical  formulae  are  contained  in  this  report. 

A.  TEST  CONFIGURATION 

Figures  22  and  23  illustrate  the  test  configurations  for  the  round  rod  and  curved  beam 
supports.  (The  configuration  for  the  squirrel  cage  supports  is  similar  to  that  of  the  round  rod 
supports.)  These  fixtures  simulate  the  mounting  of  the  supports  in  an  actual  engine  installation. 
Massive  “plugs”  were  mounted  inside  the  ID  of  the  round  rod  and  squirrel  cage  supports  when 
necessary  to  simulate  the  bearing  outer  races  and  to  prevent  ovalization  of  the  support  during 
testing.  Support  fixtures  were  bolted  to  a massive  steel  table  to  assure  end  rigidity.  Figures  24  and 
25  are  photographs  of  the  squirrel  cage  and  curved  beam  support  test  rigs,  respectively. 
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Figure  25.  Curved  Beam  Support  Test  Rig 


A force  applied  to  the  flexible  support  using  an  electromechanical  shaker  was  measured  by 
the  attached  load  cell  and  signal  amplifier  over  a frequency  range  of  10  to  2,000  Hz;  the  response 
was  measured  by  an  accelerometer  mounted  directly  opposite  the  shaker.  Plots  were  then  made 
of  mechanical  impedance  versus  frequency;  the  stiffness  and  damping  values  were  then  directly 
extracted  from  the  impedance  plots.  For  a vibrating  spring-mass-damper  system  being  excited  by 
a sinusoidal  force  at  a force  at  .<*  frequency  less  then  one-half  the  natural  frequency,  the  primary 
resistance  force  is  due  to  the  spring  force.  At  the  natural  frequency  the  resistance  force  is  due  to 
damper.  Taking  advantage  of  these  characteristics  of  a dynamic  system,  the  spring  and  damping 
rates  of  the  bearing  supports  were  determined  from  the  mechanical  impedance  plots.  These 
values  were  designated  as  the  experimental  dynamic  stiffness  and  damping  terms. 

B.  RE8ULT8 

1.  For  the  squirrel  cage  supports,  the  variation  between  the  experimental 
dynamic  stiffness  and  the  calculated  static  stiffness  ranged  from  -6%  to 

+8%. 

2.  For  the  round  rod  supports,  the  variation  between  the  experimental  dynamic 
stiffness  and  the  calculated  static  stiffness  ranged  from  -15%  to  -51%. 

3.  For  the  curved  beam  supports,  the  variation  between  the  experimental 
dynamic  stiffness  and  the  calculated  static  stiffness  ranged  from  +5%  to 
+27%. 
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4.  The  damping  from  the  squirrel  cage  and  round  rod  supports  was  negligible. 

5.  The  damping  from  the  curved  beam  supports  was  significant,  ranging  from 
4.1  to  10.0  ib-sec/in. 

Note  that  in  Results  1,  2,  and  3 the  experimental  dynamic  stiffness  was  compared  to  the 
calculated  static  stiffness.  This  is  becauce  the  static  stiffness  measured  at  low  frequencies  was 
found  to  equal  the  dynamic  stiffness. 

A summary  of  the  experimental  and  analytical  results  is  provided  in  Table  3. 

The  experimental  values  shown  for  structural  damping  are  fairly  consistent..  Yet  the 
minimum  and  maximum  damping  values  for  the  round  rod  supports  differ  by  a factor  of  ten. 
Following  is  a possible  explanation:  Inside  each  of  the  round  rod  supports  is  a sleeve  up  i which 
the  piston  rigs  are  mounted.  In  the  assembly  process,  it  is  possible  that  the  sleeve  may  hevo  been 
brared  off-center  such  that  rubbing  occurred  between  the  sleeve  and  the  inner  surface  of  the  main 
body  of  the  support.  This  would  explain  both  the  large  variation  in  structural  damping  and  the 
large  amount  of  damping  that  two  of  the  round  rod  supports  had  when  compared  with  the  squirrel 
cage  supports. 

Ncte  the  high  amount  of  damping  that  the  curved  beam  supports  exhibited.  This  was  due 
to  the  following:  As  the  beam  was  deflected  by  the  load  in  the  y direction,  the  chord  AB  changed 
in  length  (Figure  26).  As  the  chord  changed  in  length,  the  supporting  ends  of  the  beam  moved 
circumferentially  along  the  supporting  surface,  causing  frictional  damping. 

C.  SUPPORT  STIFFNESS  FORMULA 

For  round  rod  supports,  the  following  formula  should  be  used  to  predict  the  static  stiffness: 

K u 12  NE  y d4 

64L*  + 112.32d:L 


For  squirrel  cage  supports: 


R = 


NEbt 


(- 


i 


i 


L*/t*  + 13176  L*/b*  + 13L/6 


Here 


K is  the  stiffness 
b is  the  leg  width 
t is  the  leg  thickness 
L is  the  length  of  the  leg  (see  Figure  27) 

Castigliano’s  Theorem  and  Curved  Beam  Theory  yield  the  following  equations  for 
calculating  the  support  stiffness  of  a 3-beam  support: 

„ „„„  Et*h  , r,  -f  r, 

K = 0.331  — — , where  » — - — , t » r,  - rt 

r\v  2 


44 


Figure  26.  Curved  Beam  Geometry 


Squirrel  Cage 


1 

b 

T 


Round  Rod 


Figure  27.  Beam  Dimensions 


Hera  K is  the  support  stiffness,  and  E is  the  material  modulus  of  elasticity.  The  support 
dimensions  are  shown  in  Figure  28. 

NOTE:  This  formula  becomes  increasingly  inaccurate  as  r,v/t  becomes  less  than  10. 
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Figure  28.  Curved  Beam  Support  Dimensions 


D.  CONCLUSIONS 

Existing  beam  formulas  for  calculating  the  Btiffness  of  flexible  bearing  supports  have 
satisfactory  accuracy  for  estimating  support  performance.  If  a maximum  error  of  about  20% 
cannot  be  tolerated,  then  a more  precise  modeling  technique  such  as  a finite  element  method 
should  be  used. 
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SECTION  IV 


SYNCHRONOUS  WHIRL  RIQ  — ANALYSIS  VERIFICATION  (TASKS  V,  VI,  AND  VII) 


A comprehensive  program  was  conducted  to  both  analytically  predict  and  experimentally 
verify  the  dynamic  response  of  a viscous-damped  flexible  rotor  system  subjected  to  both  normal 
and  abnormal  imbalance  conditions.  The  test  model  simulated  future  front-drive  power  turbine 
configurations;  i.e.,  it  consisted  of  a long  slander  shaft  with  flexible,  viscous  damped  bearing 
supports  and  an  overhung  drive  turbine. 

This  program  was  designed  to  determine,  by  parametric  testing,  the  optimum  combi  ation 
of  stiffress  and  damping  for  subcriticai  and  supercritical  shafting.  The  optimum  design  is  defined 
as  one  which  simultaneously  limits  bearing  loads  tc  low  levels,  controls  rotor  motion  to  prevent 
shroud  contact,  and  absorbs  large,  suddenly  applied  imbalance  without  failure. 

Analytical  predictions  employed  were  made  using  the  viscous  damper  analysis  methods 
substantiated  in  Task  IQ, 

A.  EXPERIMENTAL  RIG  DESIGN 

The  synchronous  reponse  rig  (Figures  29  through  31)  consists  of  a main  shaft  with  two  disks 
mounted  near  one  end  to  simulate  * front-drive  power  turbine  rotor.  The  shaft  is  provided  with 
five  balance  planes.  This  allows  for  multiplane  balance  of  the  rig  as  well  as  provides  controlled 
imbalance.  The  shaft  measures  14.2  in  from  bearing  to  bearing  with  a diameter  of  0.9  in.  and  a 
weight  of  20.5  lb.  The  front  end  of  the  shaft  incorporates  a preloaded  duplex  bearing  and  viscous 
damper.  However,  this  damper  was  inactive  during  the  test  to  isolate  the  effects  of  the  turbine 
and  damper  being  studied.  The  aft  end  of  the  shaft  is  supported  by  a roller  bearing.  The  roller 
bearing  is  also  equipped  with  a viscous  damper.  The  damper  length  and  diametei  are  1.2  in.  and 
4.0  in.  respectively.  Damper  parameters,  such  as  oil-film  thickness,  bumper  clearance,  and 
bearing  support  spring  rate,  were  varied  during  test. 
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Figure  30.  Synchronous  Response  Rig,  Three-Quarter  Rear  View 


Simulator - 


Figure  31.  Synchronous  Response  Rig  Schematic  Diagram 


Cantilevered  from  the  roller  bearing  is  the  drive  turbine  from  which  the  blade  loos  simulator 
disk  is  mounted,  as  shown  in  Figures  31  and  32.  The  blade  loss  mechanism  consist?  of  two  release 
arms  pivoted  on  a clevis.  The  arms  were  tied  together  holding  the  1/3-os  in.  weights  which  were 
released.  The  weights  were  released  by  cutting  the  tie  string  via  a blade  attached  to  a pneumatic 
actuator  (Figure  33). 


Figure  32.  Synchronous  Response  Rig  Showing  Blade  Lo3S  Simulator  bisk 


B.  R)G  TEST 
1.  T*st  Program 

The  test  instrumentation  is  depicted  in  Figure  34.  It  consisted  of  six  pairs  of  proximity 
probes  to  determine  whirl  orbit  at  selected  shaft  locations.  In  addition,  a set  of  accelerometers 
and  proximity  probes  were  used  to  measuij  the  response  of  the  viscous  damper  itself. 

A strain  gage  located  on  the  outer  damper  housing  measured  the  load  transmitted  through 
the  damper.  The  accelerometers  placed  on  th«  calibrated  squirrel  cage  bearing  support  allowed 
for  the  force  measurement,  at  the  bearing  support.  The  resultant  force  determined  the  magnitude 
of  the  bearing  load. 

All  data  was  monitored  via  oscil'.oscope  and  recorded  on  magnetic  tape  for  subsequent  data 
reduction.  A photograph  of  the  test  stand  controls  and  instrumentation  is  shown  in  Figure  35. 

The  test  program  is  outlined  in  Figure  36.  An  initial  balance  of  the  rig  was  followed  by 
running  the  rig  using  three  predetermined  imbalances.  For  the  three  imbalances  there  were  two 
changes  in  the  viscous  damper  film  thickness  followed  by  one  change  in  bearing  support  spring 
rate  and  a change-back  to  the  initial  film  thickness.  This  series  was  followed  by  two  supercritical 
and  two  subcritical  blade  loss  simulations,  incorporating  two  different  bumper  configurations. 
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Figure  33.  Blade  Loss  Release  Mechanism 


2.  Test  Results 

The  objective  of  the  synchronous  response  rig  test  was  to  generate  parametric  experimental 
data  to  optimize  the  dynamic  response  of  the  rotor  system  under  normal  and  abusive  imbalance 
conditions  such  as  the  loss  of  a blade.  Parameters  such  as  bearing  load  and  mtor  displacement 
were  measured  as  a function  of  imbalance,  shaft  speed,  rotor-to-shroud  clearance,  bearing 
damper  support  flexibility,  and  damper  clearance.  Analytical  correlation  of  these  parameters 
with  the  experimental  data  is  also  shown. 


The  results  of  the  parametric  rig  tests,  builds  1 through  7,  are  presented  in  three  categories: 
effect  of  imbalance  and  damper  clearance,  effect  of  bearing  support  stiffness  and  damper 
clearance,  and  effect  of  bumper  clearance.  The  dynamic  responses  of  four  key  rotor  parameters 
are  presented  as  a function  of  these  three  categories.  The  key  parameters  consist  of  the  maximum 
dynamic  response  at  the  critical  speed  of  the  midspan  shaft  deflection,  turbine -shroud  relative 
deflection,  bearing  flexible  support  load,  and  damper  load. 


Blade  Loss  Disk 
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Bearing  Support 


Proximity  Probes  at  0°  Deg 

Damper  Proximity  Probes  at  0°  and  90 

Accelerometers  at  0J  and  90 
Strain  Gage  Load  Cells  at  0°  and  90 


\ Shaft  Proximity 


Probes  at  0° 
and  90° 


Figure  34.  Synchronous  Rig  Instrumentation 


Figure  35.  Rig  Test  Stand  Controls  and  Instrumentation 
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Figure  36.  Task  VU  Test  Program 
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The  effects  of  imbalance  and  damper  clearance  on  rotor  dynamic  response  are  shown  in 
Figure  37,  which  summarizes  the  results  of  rig  builds  1,  2,  and  3.  As  expected,  the  rotor-bearing 
response  increased  with  increased  imbalance;  however,  the  response  was  minimized  with  the 
largest  damper  clearance.  This  indicates  that  the  higher  damper  forces  generated  with  the  lower 
clearances  tend  to  limit  the  orbital  motion  in  the  damper,  thus  forcing  more  strain  energy  in  the 
whirling  shaft  and  reducing  the  effective  damping.  Although  both  the  bearing  loading  loads  and 
midshaft  deflections  increase  with  decreased  clearance,  the  turbine  shroud  deflection  was 
reduced  slightly  with  the  0.006-in.  clearance  damper  because  the  turbine  is  located  near  the 
damper  where  the  orbital  motion  is  forced  to  be  limited.  However,  the  0.010-in.  clearance  damper 
reduced  dynamic  response  of  all  four  key  rotor  parameters.  It  should  be  noted  that  as  the 
imbalance  increases,  the  damper  load  increases  at  a faster  rate  than  the  hairpin  load.  This  is  an 
indication  of  the  characteristic  nonlinear  force  effects  of  the  squeeze  film  damper.  As  the  bearing 
deflection  amplitude  increases,  the  damper  transmits  more  and  more  of  the  bearing  load, 
reducing  the  load  transmitted  through  the  hairpin. 

The  effects  of  bearing  support  (hairpin)  stiffness  and  damper  clearance  on  rotor  dynamic 
response  are  shown  in  Figure  38,  which  summarizes  the  results  of  rig  builds  2,  3,  4,  and  5.  The 
most  obvious  result  is  that  rotor  response  is  reduced  by  reducing  the  bearing  support  (hairpin) 
stiffness.  Reducing  the  stiffness  allows  an  increase  in  orbital  motion  in  the  squeeze  film  damper, 
which  increases  effective  damping  and  reduces  vibration.  The  increase  in  damper  load  with 
decreased  hairpin  stiffness  also  indicates  the  increased  damper  motion.  Again  the  minimum  rotor 
response  is  observed  when  the  largest  damper  clearance  (0.010  in.)  is  used.  These  data  clearly 
indicate  that  soft  bearing  supports  in  conjunction  with  wide  damper  clearances  are  a good 
combination  for  this  flexible  rotor  rig. 

The  effects  of  bearing  bumper  clearance  on  the  dynamics  of  a flexible  rotor  under  transient 
blade  loss  conditions  for  both  subcritical  and  supercritical  speeds  are  sho  vn  in  Figure  39.  These 
plots  are  slightly  different  from  the  previous  ones  in  that  the  peak  transient  response  is  displayed, 
as  opposed  to  the  peak  steady -state  response  on  the  previous  plots.  The  most  significant  result  of 
the  blade  loss  simulations  is  that  the  peak  rotor  transient  response  is  reduced  with  increased 
bumper  clearance.  This  indicates  that  as  the  bearing  support  strikes  the  bumper,  high  loads  are 
transferred  to  the  rotor  and  at  this  point  the  effectiveness  of  the  damper  is  truncated  because  no 
further  excursion  into  the  damper  film  is  permitted  by  the  bumper.  When  the  larger  bumper 
clearance  is  used,  more  effective  squeeze  film  damping  is  achieved  and  the  bumper  contact  loads 
are  reduced.  This  is  also  indicated  by  the  decreasing  hairpin  load  with  decreased  bumper 
clearance,  because  once  the  bumper  is  struck,  the  hairpin  is  relieved  of  an  additional  load  and 
further  load  is  transmitted  through  the  damper  support.  Note  the  increased  damper  support  load 
with  decreased  bumper  clearance.  Another  observation  of  the  blade-loss  test  results  is  that  the 
rotor  response  amplitudes  are  always  higher  for  supercritical  rotor  operation  than  for  subcritical 
operation. 

The  scope  of  the  test  program  was  slightly  expanded  in  order  to  investigate  the  effects  of 
decelerating  a supercritical  rotor  through  its  critical  speed  after  a blade  loss.  The  test  was 
performed  at  the  request  of  the  Army  Program  Manager  to  verify  that  no  severe  damage  would 
occur  to  an  engine  with  supercritical  speed  shafting  under  blade  loss  conditions.  This  test  differed 
from  the  previous  supercritical  blade  lose  test  in  that  the  rotor  was  not  rebalanced  by  jettisoning 
a second  balance  weight  before  deceleration.  Figure  40  shows  the  midshaft  vibration  amplitude 
of  the  rotor  acceleration  through  the  critical  speed  with  only  residual  imbalance,  then  sudden 
blade  loss  and  the  deceleration  back  through  the  critical  speed  with  the  blade  loss  imbalance. 
Although  the  rotor  rubbed  bearing  compartment  seals,  no  structural  damage  to  the  rig  was 
incurred. 
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Figure  40.  Midshaft  Vibration  Amplitude  vs  Speed 

In  summary,  the  parametric  test  program  has  been  successful  in  giving  a much  better 
understanding  of  how  variations  in  the  most  significant  of  rotor  design  parameters  affect  rotor 
response  to  synchronous  dynamic  excitation.  The  main  results  were  that  increased  damper  and 
bumper  clearance,  and  decreased  support  (hairpin)  stiffness,  minimized  the  dyne-  espouse  of 
the  flexible  shaft  rig.  In  engine  applications,  however,  tnings  such  as  maneuve:  umpressor 

surge  load  deflections  limit  how  “soft”  the  bearing  support  can  be  and  how  la.  the  bumper 
clearance  can  be.  The  best  combination  of  these  parameters  for  optimum  rotor  performance  is  a 
function  of  the  particular  engine  durability  and  efficiency  requirements.  Accurate  analytical 
models  of  the  rotor  system  are  the  most  viable  means  of  determining  the  best  trade-off  of  all 
variables.  The  next  section  contains  the  comparison  of  analytical  predictions  of  commonly  used 
rotor  dynamics  commuter  programs  with  the  flexible  rotor  rig  test  data. 


C.  DATA  CORRELATION  WITH  ANALYTICAL  MODEL 


Two  basic  types  of  rotor  dynamic  response  models  we  re  used  to  predict  the  dynamic  response 
of  the  synchronous  whirl  rig.  The  first  model  was  used  to  predict  the  rotor  dynamic  response  to 
normal  steady-state  imbalance  loads  as  were  used  in  builds  1 through  6 of  the  test  program.  This 
model  is  the  commonly  uned  transfer  matrix  forced  response  computer  program  similar  to  that  of 
Reference  10.  The  second  model  was  used  to  predict  the  transient  blade  loss  dynamics  cf  the 
synchronous  rig,  builds  6 and  7 of  this  test  program.  This  computer  program  used  the  v/ell-known 
modal  technique  for  calculating  the  transient  response  of  large  dynamic  systems  and  is  similar  to 
that  of  Reference  11.  Each  of  these  mathematical  models  used  the  squeeze  film  damper 
mathematical  model  developed  in  Tank  I of  this  contract  to  predict  the  damper  forces. 

10  Lund,  J.  W.,  and  Or^utt,  F.  K.,  "Calculations  and  Experiment*  on  the  Unbeltm-e  Response  of  a Flexible  Kotor." 
ASME  Paper  No.  67-Vibr-27. 

11  Dennis,  A.  -J.,  Erikseoi.  R.  H.,  and  Seitelman,  L.  H.,  "Transient  Response  Analysis  of  Pumped  Rotor  Systems  by  the 
Normal  Mode  Mefhod,"  ASME  Paper  No.  75-GT-58. 
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The  first  step  in  predicting  the  dynamic  response  of  a rotor  system  is  to  calculate  the  critical 
speed.  The  critical  speed  prediction  aided  in  determining  the  severity  of  a natural  rotor  mode  and 
the  particular  speed  to  concentrate  on  in  predicting  the  steady-state  forced  response  of  the  rotor. 
Calculation  of  tne  critical  speed  is  also  necessary  because  the  predicted  mode  shapes  and  natural 
frequencies  are  used  by  the  transient  response  model  to  predict  the  rotor  deflections  and  forces 
due  to  transient  loadings  such  as  blade  loss.  Figure  41  shows  the  cross  section  of  synchronous 
response  rig  and  the  corresponding  beam  model  which  was  used  to  predict  the  critical  speeds  of 
the  rig.  The  Hist  three  rig  critical  speeds  are  shown  in  Figure  42  along  with  the  corresponding 
relative  mode  shapes  and  rctor  strain  energies.  The  5000-rpm  critical  speed  is  the  rotor  mode  that 
is  being  investigated  throughout  this  test  program.  The  91%  rotor  strain  energy  in  this  mode 
indicates  a very  severe  critical  speed,  because  very  little  effective  damping  can  be  produced  by 
the  squeeze  film  damper  when  so  little  of  the  vibratory  energy  exists  in  the  support  system. 
Compared  to  current  small  engine  standards,  this  program  is  an  aggressive  step  forward  in 
supercritical  speed  rotor  development. 

The  next  step  in  predicting  the  dynamic  performance  of  a rotor  system  is  to  calculate  the 
system  response  amplitude  vs  speed  as  a function  of  the  applied  load.  The  applied  load  for  builds 
1 through  5 of  the  synchronous  whirl  rig  test  program  was  an  imbalance  placed  at  the  blade  loss 
disk.  A representative  comparison  of  the  predicted  rig  response  to  the  measured  response  is  shown 
in  Figures  43  through  46.  Shaft  vibration  sensitivity  and  bearing  load  sensitivity  were  chosen 
because  of  their  influence  on  engine  performance  degradation  and  durability.  The  speed  at  which 
the  peak  amplitude  of  response  occurred  was  predicted  with  good  accuracy,  and  the  average 
deviation  of  measured  and  predicted  peak  amplitudes  was  about  17%.  Considering  a combination 
of  measurement  and  computation  error,  steady-state  test  data  are  in  good  agreement  with 
prediction. 

Computation  error  is  calculated  in  terms  of  system  energy  error.  This  error  is  a function  of  the 
number  of  significant  digits  that  can  be  maintained  in  a computer,  and  this  numerical  error  can 
accumulate  in  some  programs  to  the  point  where  the  predictions  have  considerable  deviation 
from  the  exact  values.  The  transfer  matrix  numerical  method  used  in  this  report  and  in  most 
steady-state  rotor  dynamic  response  programs  is  susceptible  to  accumulative  round-off  error.  To 
assess  the  amount  of  numerical  error  that  has  been  accumulated,  comparison  is  made  between 
the  energy  input  to  the  dynamic  system  through  the  imbalance  load  and  the  energy  dissipated  out 
through  the  dampers.  Ideally,  these  two  quantities  are  equal.  The  percent  difference  between  the 
energy  “in”  and  the  energy  “out”  is  a good  measure  of  the  numerical  error  accumulated  in 
predicting  the  system  dynamic  response.  An  average  energy  error  of  1-3%  was  observed  for  the 
synchronous  whirl  rig  imbalance  response  predictions.  Energy  errors  of  this  magnitude  are 
usually  incurred  when  predicting  the  response  of  a very  flexible  rotor  such  as  the  one  being 
studied,  and  the  resulting  prediction  errors  can  be  of  significant  magnitude.  Therefore,  numerical 
methods  that  eliminate  the  accumulation  of  round-off  error  should  be  applied  to  flexible  rotor 
dynamic  predictions. 
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Figure  41.  Synchronous  Whirl  Rig  Model 
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Figure  45.  Shaft.  Sensitivity  With  Soft  Bering  Support  (Build  4) 
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Figure  46.  Bearing  Load  Sensitivity  With  Soft  Bearing  Support  (Build  4) 
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The  final  two  tests  cf  this  seven-part  test  program  covered  transient  blade  loss  simulation  at 
both  subcritical  and  supercritical  speeds.  The  measured  data  and  predicted  response  of  the 
build  6 run  with  the  large  bumper  clearance  are  shown  in  Figures  47  and  48,  the  first  figure 
showing  the  subcritical  blade  loss  results  and  the  latter  figure  the  supercritical  remits.  Figures  49 
and  50  show  the  corresponding  blade  loss  response  plots  with  a small  bumper  clearance.  All  the 
responses  are  plotted  vs  time  from  the  point  of  blade  loss  to  steady  state.  The  measured  data  are 
shown  as  a complete  time  trace,  and  the  predicted  results  are  shown  as  an  envelope  of  peaks  and 
are  superimposed  on  the  measured  data.  Note  that  both  the  measured  and  predicted  curves  show 
the  beat  frequency  of  excitation  beating  with  the  natural  frequency  of  the  rotor.  The  measured 
and  predicted  curves  are  similar  in  character,  but  in  all  cases  the  measured  response  is  slightly 
higher  than  the  predicted.  Here  again,  normal  measurement  error  and  computational  error  can 
easily  combine  to  give  the  differences  shown. 
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Figure  47.  Subcritical  Blade  Loan  Response  (Lorg^i  Bumper  Clearance) 
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SECTION  V 


NON8YNCHRONOUS  WHIRL  RIQ  — ANALY8I8 
FORMULATION  ANO  VERIFICATION  (TASKS  VIII  AND  IX) 


A.  ANALYTICAL  MODEL 
1.  Description 

Spline  couplings  are  commonly  used  in  high  speed  rotating  machinery  to  transfer  power 
between  adjacent  shafts  or  within  whet  is  nominally  one  shaft  for  ease  of  assembly/disassembly. 
Normally,  the  two  halves  of  these  couplings  are  rigidly  aligned  and  fixed  with  respect  to  each 
other  by  means  of  a pilot  so  that  no  relative  motion  is  possible  between  adjacent  spline  teeth 
except  that  due  to  elastic  deflection  resulting  from  transmitted  torque.  This  is  not  necessarily  the 
case  for  unpiioted  spline  couplings.  For  this  class  of  couplings,  alignment  in  both  radial  and 
angular  directions  is  provided  solely  by  coupling  geometry  and  close  tolerances  on  the  mating 
parts.  In  addition,  no  axial  restraint  exists.  A particular  application  of  an  unpiloted  spline 
coupling  might  be  the  connection  between  the  power  turbine  shaft  and  transmission  input  shaft 
in  a helicopter  power  train.  For  this  case,  the  female  coupling  half  is  supported  on  bearings  while 
the  male  half  is  integral  with  the  turbine  shaft  which  is  only  supported  on  bearings  at  its  far  end. 
Since  the  coupling  is  not  solidly  locked  together,  relative  motion  along  several  axes  may  occur 
between  coupling  helves.  The  forces  generated  by  the  friction  between  moving  surfaces  give  rise 
to  forces  within  the  coupling  which  tend  to  destabilize  the  rotor  system  and  can  indeed  cause 
instability  if  certain  “favorable"  conditions  (for  example,  insufficient  external  damping)  exist. 

This  problem  has  occurred  in  the  past  and  is  the  subject  of  a paper  by  Williams  and  Trent 
(Reference  12),  who  explored  the  effects  of  asymmetry  and  nonlinearity  in  the  supporting 
structure  on  spline-coupling-induced  instability.  Their  study,  however,  assumed  a rigid  rotor  and 
used  an  analog  computer  to  simulate  the  rotor  bearing  system. 

This  report  documents  an  analysis  cf  unpiloted  spline  couplings.  The  coupling  model  thuo 
derived  can  be  incorporated  into  a rotor  dynamics  computer  program  to  determine  the  stability 
of  a rotor-bearing-spline  coupling  system. 

In  this  report,  the  spline  coupling  is  modeled  as  a section  of  the  shaft  system  having  internal 
damping.  (For  a further  discussion  of  internal  damping  and  its  effects  on  stability,  see  Reference 
13.)  This  model  may  then  be  incorporated  as  one  portion  of  a rotor-bearing  system  model.  The 
system  stability  may  then  be  studied  by  computing  the  complex  eigenvalues  or  damped  natural 
Kjuencies  of  the  system.  Each  damped  natural  frequency  is  a complex  number  wherein  the 
imaginary  part  represents  a frequency  of  oscillation  of  the  system  and  the  real  part  indicates 
stability.  Positive  or  negative  values  for  a real  part  indicate  that  the  system  will  be  unstable  or 
stable,  respectively,  in  that  mode.  A stability  map  (plot  of  real  part  against  some  parameter  such 
as  peed)  or  a Nyquist  plot  (real  against  imaginary  part  for  various  conditions)  can  then  be  drawn 
up  covering  the  modes  and  operating  conditions  of  interest.  An  overall  picture  of  the  stability  of 
the  system  can  thus  be  generated. 


12  Williams,  R.,  and  Trent,  R„  "The  Effect  of  Non-Linear  Asymmetric  Supports  on  Turbine  Engine  Rotor  Stability,” 
SAE  700320. 

13  Lund,  J.  W.,  "Stability  and  Damped  Critical  Speed*  of  ■ Flexible  Rotor  in  Fluid  Film  Bearings,"  Journal  of 
Engineering  for  Industry,  Trans.  ASME,  Series  B,  Vol.  96,  No.  2,  May  1974,  pp.  609-617. 
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For  each  point  in  the  analysis,  the  spline  coupling  model  requires  four  coefficients:  two 
stiffness  and  two  damping  terms.  These  coefficients  are  derived  in  Sections  V.A.2.a  aid  V.A.2.b, 
respectively,  and  are  functions  of  spline  geometry  and  coefficient  of  friction.  Section  V.A.2.C 
suggests  an  approach  for  stability  calculation  and  offers  guidance  on  the  selection  of  coefficient 
of  friction. 

2.  Method  of  Solution 
a.  Spllna  Coupling  Stiffness 

The  stiffness  of  an  unpiloted  spline  coupling  may  be  calculated  by  first  rinding  the  stiffness 
of  the  individual  coupling  teeth  for  various  loading  conditions.  The  overall  coupling  stiffnesses  in 
the  lateral  and  angular  directions  are  then  found  by  properly  combining  the  individual  tooth 
stiffnesses.  For  this  analysis,  the  following  assumptions  will  be  made: 

• Except  for  a small  localized  area  about  the  base  of  each  tooth,  the  material 
on  which  the  teeth  are  mounted  is  rigid.  This  implies  that  the  deflection  of 
the  coupling  results  entirely  from  deflection  of  the  teeth  in  rigid  supports. 

• All  deflections  are  very  small  and  all  deflecting  parts  behave  elastically. 

• Where  a given  stiffness  is  nonlinear  (nonproportional  to  deflection),  it  will  be 
linearized  about  its  value  under  steady-state  operating  conditions. 

• Errors  in  tooth  location  and  profile  are  small.  Each  tooth  is  in  line  ccntact 
with  its  mate,  and  all  carry  about  the  same  load. 

• The  forces  on  the  coupling  teeth  are  due  predominantly  to  the  transmitted 
torque  and  forces  arising  from  the  misalignment  itself  (such  as  elastic 
restoring  forces).  Other  forces,  such  as  those  resulting  from  friction-induced 
moments,  are  small. 

The  analysis  for  the  individual  tooth  stiffness  is  similar  to  that  in  Reference  14,  except  for 
the  Hertzian  contact  deflection  term,  which  may  be  found  in  Reference  15. 

Spllna  Tooth  Stiffness  — Tangential  Direction 

The  stiffness  of  a spline  tooth  may  be  found  by  computing  the  deflection  8 under  load  of  a 
single  tooth  as  the  result  of  each  of  four  mechanisms: 

da  = displacement  due  to  bending  as  a cantilever  beam 

8b  - displacement  due  to  shear  deformation  as  a cantilever  beam 

8r  = displacement  due  to  tooth  rotation  as  a rigid  body  in  the  supporting 
structure  at  its  base 

~ displacement  due  to  Hertzian  contact  deformation  of  the  tooth  surface 


14  Laskirs,  1.,  et  al„  “Analysis  os'  Noise  Cienerated  by  UH-1  Helicopter  Transmission,"  USAAVLABS  Technical  Report 
68-41,  June  1968,  Available  through  Clearinghouse,  Springfield,  Va.  22151,  Report  No.  AD  676468. 

15  Walowit,  J,  A.,  and  Anno,  J.  N„  “Modem  Developments  in  Lubrication  Mechanics.’  Applied  Science  Publishers. 
Ltd.  London  1976. 
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Outplacement  due  to  transverse  shear  of  the  supporting  structure  and  radial  compression  of 
the  tooth  are  neglected  (the  latter  for  tangential  stifftiess  only).  The  displacements  will  then  be 
used  to  find  tooth  compliance  resulting  from  each  mechanism.  The  individual  compliances  are 
then  simply  added  and  the  result  is  inverted  to  yield  the  tooth  stiffness. 

To  calculate  SB  and  i8.  the  spline  tooth  will  be  treated  as  a sequence  of  transverse  sections, 
each  of  uniform  rectangular  cross  section.  The  depth  of  each  segment  is  the  average  cf  the  depths 
at  the  ends  of  the  segment,  as  shown  in  Figure  51.  In  this  figure,  yk  represents  the  distance  from 
the  tooth  surface  to  its  centerline  at  radius  rk,  rp  is  the  pitch  radius,  P is  the  applied  load,  and  * 
is  the  angle  from  the  horizontal  of  the  applied  load. 


Deflection  Due  to  Bending 

The  second  moment  of  area  of  a tooth  cross  section  at  isdius  rk  is  expressed  as: 

h = -jj  «(2  yk)*  = J- 1 y{  (28) 

where!  is  the  length  of  the  tooth  face.  At  radius  rk.p  at  the  other  end  of  the  tooth  segment,  the 
second  moment  of  area  is 

ii.,-  j1  y e-  (29) 

For  the  assumed  cross  section,  the  mean  of  the  two  second  moments  of  area  will  be  used: 
l = — -2-~—  (30a) 

or 

l = j (yt  + yg_k)  (30b) 
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The  average  a.ea  A-  of  the  segment  can  be  found  in  a similar  manner,  resulting:  in 

A#  - t (y*  + y»-,)  (31) 

The  height  of  segment  k may  be  found  from  the  radius  at  cither  end  of  the  segment: 

hte  « r*  - r*_.  (32) 

Each  of  the  above  equations,  (30b)  (31),  and  (32),  is  applied  to  every  segment  on  the  tooth, 
both  on  the  involute  curve  and  at  th?  tooth  root. 

The  distance  of  the  k*h  segment,  at  its  cuter  end,  from  the  point  of  load  application  is  fouod 

to  be 


S*  = ?p  - rk  (33) 

where  the  face  width  of  the  spline  is  large  compared  with  the  tooth  thickness:  the  extra 
width  contributes  an  added  stiffening  effect  to  the  teeth.  This  may  be  incorporated  into 
the  bending'  deflection  equation  tkrougn  use  of  the  factor  (l-rf),  where  » is  Poisson  a ra'.io. 
Normally,  this  appears  in  the  numerator  while  the  modulus  of  elasticity  E appears  in  the 
denominator.  For  convenience,  the  (I-**)  factor  will  be  applied  as  a correction  to  the 
modulus  E to  yield  a corrected  value  E*,,  which  is  defined  as  follows: 

E,  - (341 

Should  the  spline  teeth  be  relatively  narrow  compared  with  their  thickness,  the  correction  factor 
should  be  dropped. 

In  Figure  51,  the  lead,  P,  acts  in  two  ways  to  load  the  cantilevered  spline  tooth: 

1.  A concentrated  load  equal  to  the  transverse  component  of  the  applied  load, 

P cos  ♦. 

2.  A moment  load  equal  to  the  radial  component  of  the  applied  load  times  one- 
half  the  thickness  of  the  tooth  at  the  point  of  load  application,  Pyp  sin  ❖. 

In  tooth  bending,  each  of  these  load  conditions  causes  a transverse  displacement  and 
rotation  of  the  tooth  centerline  under  the  point  of  load  application.  These  two  displacements  will 
bo  evaluated  for  each  of  the  cwo  load  conditions  and  then  combined. 

Figure  52  portrays,  in  its  simplest  form,  a tooth  segment  subjected  to  a transverse  load 
applied  some  distance  away.  The  rest  of  the  tooth  serves  as  a rigid  support  or  a rigid  extension. 
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Tooth  Segment 


Figure  52.  Bending  of  a Tooth  Segment  Subjected  to  a Transverse  Load 


From  standard  beam  deflection  analysis,  the  transverse  deflection  j,  and  centerline  rotation 
>1  at  the  load  are 
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The  two  equations  can  be  rewritten  using  the  terminology  of  the  main  analysis,  so  that  for  a 
segment  k, 


(5Bk)l  - 


Ply  coe'k 
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(h{  + 3 S*  hk 
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(7bJi  = 
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(38) 


The  other  case  in  tooth  bending,  shown  in  Figure  53,  considers  a tooth  segment  subjected  to 
a moment  load  applied  some  distance  away,  with  the  rest  of  the  tooth  serving  as  a rigid  support 
or  extension. 


Tooth  Segment 


Figure  53.  Bending  of  a Tooth  Segment  Subjected  to  a Movement 
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Again,  from  beam  deflection  analysis,  the  transverse  deflection  6,  and  centerline  rotation 
7,  at  the  load  are 


Mh 
2 El 


(h  + 2S) 


(39) 


7i  = “ 


Mh 

El 


(40) 


The  minus  sign  is  applied  because,  with  a moment  applied  in  the  direction  Bhown,  the  deflections 
are  opposite  those  described  in  the  previous  case.  Again,  Equations  (39)  and  (40)  may  be 
rewritten  using  the  main  terminology  and  thus  become 


(aBk)s  ~ 

(7Bk)f  = 


-Pyp  hk  sin* 
2E„  Ik 

-Pyp  hk  cos* 

E„Ik 


(hk  + 2 Sk) 


(41) 

(42) 


The  deflections  for  each  tooth  segment  are  next  summed  by  rewriting  Equations  (37),  (38), 
(41),  and  (42)  with  the  summation  covering  all  segments  from  the  ^ase  of  the  tooth  to  the  point 
of  the  applied  load. 
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These  displacements  are  combined  by  considering  each  load  condition  separately.  The  results 
will  then  yield  the  total  displacement  in  bending.  Figure  54  shows  the  geometry  from  which  the 
relationships  between  the  various  displacements  may  be  found. 

From  Figure  54  we  see  that 

JD  = JG  - DG  (47) 


Hence, 


JD  - 


JJ 


cos*  cos* 


(48) 


when  ie  the  displacement  of  the  tooth  outline  in  the  direction  of  the  applied  force.  Next  we  see 
that 

JG  = CC ' = &y  (49) 

where  by  ia  the  transverse  displacement  described  above.  Vector  DG  can  be  found  by  the  relations 


DG  = GG'  tan  * (60) 

GG  » GC  (ZGC'G')  (51) 

GC  = JC'  = yp  (52) 

where  yp  is  half  the  tooth  thickness  at  the  pitch  radius  and  angular  rotation  (ZGC'G')  is  found 
to  be 

y = (ZGC'G')  (53) 

which  is  described  by  Equations  (44)  and  (46).  Combination  of  Equations  (47)  through  (53)  yields 

~~  =Sr-yPy  tan  * (54) 

where  dB  is  the  total  motion  of  point  J along  the  line  of  force  due  to  bending  of  the  tooth.  Equation 
(54)  can  be  simplified  by  multiplying  through  with  cos  p. 


bB  = by  cos  ~ y„  y sin  'if 
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Figure  54.  Tooth  Pending  Geometry 
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Equation  (56)  may  now  be  used  to  combine  the  individual  displace  mounts.  For  the 
concentrated  load  condition,  Equations  (43)  (44),  and  (5t>)  are  combined  to  gei, 

. , P T coa*  ❖ v hk  „ . „ _ , „ 

(5»)i  = [ — g — 2^  -y — K 3 S„  h„  + 3 60 


coe  ❖ tin  ♦ 


yv 


K 
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(hk  + 2 S,,) 


(56) 


Similarly,  for  the  moment  load  condition,  Equations  (*'o),  (46),  and  (47)  are  combined  to 

yield 


(57) 


Lastly,  the  two  can  be  combined  co  give  the  total  deflection  due  to  bending: 

5b  = (5b)i  4-  (5j))i  \58) 

Deflection  Due  to  Shear 

The  shear  deflection  of  the  spline  tooth  as  a cantilever  beam  is  caused  by  the  same 
transverse  component  of  the  applied  load,  namely  P cos  4'.  Deflection  of  this  type  simply  offsets 
the  centerline  without  rotating  any  transverse  section.  The  total  deflection  then  becomes  the  sum 
of  the  deflections  for  each  of  the  pre  viously  defined  tooth  segments.  Figure  55  epresents  a flexible 
segment  subjected  to  shear,  with  the  remainder  of  the  tooth  serving  as  a rigid  support  or 
extension. 
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Using  standard  analytical  methods,  the  transverse  deflection  for  a segment  with  a 
rectangular  cross  section  is 

, 1*2  Ph  /(-/yi 

5 “ “ga^  <69) 

where  h is  the  segment  lengtn,  G is  the  shear  modulus,  A is  the  area,  and  1.2  is  a shear  form  factor 
for  a rectangular  cross  section.  As  noted  previously,  the  centerline  rotation  at  the  load  is 

7 = 0 (60) 

When  the  above  equations  are  rewritten  using  the  proper  terminology,  the  result  k 


1.2  Plu  coe  'if 
G A* 


I 


(yX  = 0 (62) 

Equation  (oi)  is  then  summed  over  all  segments  k to  find  the  deflection  for  the  entire  tooth 
in  shear: 

a ) - 12  p ”g  ♦ V A_  ,«,* 

(«.)-  Q 4 Ak  (63) 

To  obtain  the  displacement  in  the  direction  of  the  load,  Equation  i63)  is  substituted  into  (64): 
, 1.2  P coe*  * v ht 

«.  - c \ T>  M> 

Deflection  Due  to  Rotation  In  Elastic  Support 

Rotation  of  a spline  tooth  as  a rigid  body  in  its  elastic  supposing  structure  is  a third 
contributor  to  the  overall  deflection.  The  moment  M causing  the  fetation  results  from 
multiplication  of  the  transverse  component  of  the  applied  load  by  the  moment  arm  3/ 

M - P S,  coe  * (e5) 

where  SP  is  the  distance  from  the  applied  load  to  the  base  of  the  tooth.  Figure  56  shows  the 
loading  and  its  resultant  deflection  in  schematic  form 

According  to  Reference  16,  the  equation  for  rotation  at  the  load  is 

1.327  M 

7 ' T «rr 

which  means  that  the  displacement  at  the  load  is 


6 * y S * 


1.327  M 3 
E yi» 


18  O’Dwwl,  W.  J.,  "btnwi  and  Deflection  1b  Built-In  B**nu,”  ASME  Peper  No.  62-«VA-ie,  1W2. 
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Replacing  the  variables  in  Equation  (67)  with  those  of  the  main  analysis,  we  get 


<«*> 


1.327  P (Sp  cos  'it  - yp  sin  fr)  Sp 
E„  y>J  t 


(68) 


where  yR  is  the  half-v/idth  of  the  tooth  at  its  root  and  M has  become 

M = P (SP  cos  ♦ - yD  sin  9)  (69) 

Likewise  for  the  rotation: 


, , 1.327  P (Sp  cos  * - yp  sin  *) 

M. (.0) 

The  resulting  displacement  of  the  tooth  outline  in  the  direction  of  the  load  is  found  by 
substituting  Equations  (68)  and  (70)  into  (56): 

s = 1.327  (S„  cce  » - y„  sin  *)» 

* yi  * 


Figure  56.  Tocth  Displacer  lent  Due  to  Rotation  at  Base 


The  above  equations  have  been  derived  with  the  external  spline  tooth  in  mind.  However,  the 
same  equations  apply  for  the  internal  spline  tooth  as  long  ns  the  values  for  h and  S are  kept 
Dositive  and  as  long  as  the  summation  of  the  k segments  proceeds  between  the  base  of  the  tooth 
and  the  lead  point. 

Conversion  to  Compliance 

Equations  (58),  (64),  and  (71)  may  be  rewritten  to  express  a compliance,  or  deflection  per 
unit  loud,  which  is  defined  ns 

Q = (72) 
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Thus,  for  the  bending  case,  Equation  (56)  becomes 


Qb  = 


_1_ 

Er< 


(hi  + 3 S,  h*  + 3 S{) 


- yp  cos^  sin^ 


L 


K 

l 


(hk  + 2 Sk)  + y‘  sin1* 


(73) 


For  shear,  Equation  (64)  becomes 


Q.  = 


1.2  cos** 
G 


L 


(74) 


and  for  rotation  in  the  support,  Equation  (71)  becomes 


n - 1-327  (Sp  cosfr  - yp 

'°‘R  K v*  > 


sin'!')1 


(75) 


E,  yi 

The  total  compliance  of  the  tcoth  acting  as  e beam  is  simply  the  sum  of  the  above  compliances: 
Qbmid  = Qn  + Q,  + Qa  (76) 


The  above  holds  true  for  the  total  compliance  of  the  mating  tooth,  which  may  be  found 
independently. 

Hertzian  Contact  Deformation 

The  last  displacement  to  be  considered  in  the  tangential  direction  will  be  that  due  to 
contact  deformation.  If  the  surface  radii  of  the  mating  teeth  are  different  enough  from  each  other, 
such  that  the  width  of  the  contact  zone  is  much  smaller  then  the  depth  of  the  tooth,  the  contact 
displacement  can  be  considered  Hertzian.  If  the  radii  are  close  or  the  same,  such  that  after  initial 
“run  in”  contact  is  expected  to  occui  over  a wide  area  of  the  tooth,  then  the  deformation  is  no 
longer  Hertzian  and  will  be  treated  after  the  Hertzian  case.  For  cases  where  contact  deformation 
is  due  to  Hertzian  contact  deformation,  the  treatment  used  in  this  analysis  is  based  on  Reference 
15  for  cylindrical  surfaces  in  “line  contact.”  Figure  57  shows  two  infinitely  long  cylindrical 
surfaces  in  the  unloaded  and  loaded  conditions. 


The  two  surfaces  are  assumed  to  have  radii  of  curvature  R,  and  R*.  Unfortunately,  the 
solution  for  the  absolute  deformation  of  two  infinitely  long  cylinders  in  line  contact  is  unstable. 
However,  s solution  does  exist  for  the  change  in  length  from  the  surface  to  a point  in  each  cylinder 
remote  from  the  surface.  If  the  distances  from  the  undeformed  surfaces  of  cylinders  1 and  2 to 
remote  interior  points  are  d,  and  d,,  the  contractions  of  d,  and  d,  as  load  is  applied  are  defined 
to  be  6 m and  (Figure  57),  which  must  be  much  smaller  than  d,  and  d„  respectively. 
Additionally,  d,  and  d,  must  be  much  larger  than  the  contact  patch  width  b.  If  P is  defined  as  the 
load  per  unit  of  axial  length,  the  contraction  AH,  of  each  surface  is  given  by 


«H1  = 


2P  1-rl 

t — e , 


J- — Hi — ) 
2 (l-*,)  / 


(77) 


79 


where  the  contact  patch  width  b is  given  by 


b = 2 


(78) 


and  the  equivalent  radius  R is 


R = 


R,  R, 

Ri  + Ri 


Substitution  of  Equation  (64)  into  (78)  simplifies  it  to 


b - 2 


(79) 


(80) 


Equation  (77)  may  be  converted  to  a compliance  by  dividing  both  sides  by  the  force  P and 
substituting  Equation  (34): 


[?]-**) 


(81) 


The  above  equation  shows  that  the  compliance  is  not  constant  but  is  somewhat  dependent  on  the 
applied  load  and  the  location  from  which  the  contraction  5,  is  measured.  It  will  be  assumed  for 
this  analysis  that  distance  d,  can  be  set  equal  to  the  half- width  ypl  of  the  tooth  at  the  pitch  line. 
Since  the  dependence  of  compliance  on  load  is  not  large,  an  average  tooth  load  will  be  found  and 
then  assumed  to  be  constant  as  long  as  changes  in  displacement  are  not  too  great . 


For  cases  where  the  deformation  cannot  be  considered  Hcitzian,  the  deflection  due  to 
contact  pressure  will  be  approximated  by  testing  the  footh  as  a member  in  compression.  The 
clascicr  l equation  describing  the  deflections  of  a rod  in  compression  is 


8 


Pt_ 

AE 


(82) 


where  P is  the  applied  force,  < is  the  length  of  the  rod,  ar.d  A and  E ere  the  area  and  the  elastic 
modulus,  respectively.  The  compliance  Q is  then 


Q = 


± 

p 


t 

AE 


(83) 


For  our  case,  t will  be  taken  to  be  the  half  depth  yp  of  the  tooth  at  the  pitch  diameter,  the 
modulus  becomes  E,to  take  account  of  plane  strain,  and  the  area  is  approximated  by  length?  of 
the  tooth  face  times  the  width  wc  of  the  contact  patch.  Hence,  for  one  tooth,  i of  a pair,  L 


Qhi 


wc?E„ 


(84) 


The  subscript  for  Q will  be  the  same  as  for  the  Hertzian  case,  as  the  two  cases  are  mutually 
exclusive. 


In  summary,  the  elastic  compliance  of  a pair  of  mating  teeth  due  to  a load  acting  through 
the  tooth  flanks  is  the  sum  of  the  individual  compliances: 

Qt  = Qbi  + Qbi  + Qri  + Qhi  + Qbi  + Qhi  Qhi  + Qhi  (85) 

The  stiffness  of  the  pair  of  teeth  is  then  the  inverse  of  the  compliance: 

kT=  -J-  (86) 


Splint  Tooth  Stiffness  Radial  Direction 

The  deflection  of  a spline  tooth  in  the  radially  invard  direction  may  be  found  by  again 
splitting  the  tooth  into  segments,  finding  the  deflection  of  each  segment  to  a given  force  and 
summing  them  over  the  height  of  the  tooth.  Figure  58  shows  the  geometry  of  a spline  coupling 
tooth  loaded  radially  in  compression. 


The  equation  for  the  deflection  of  the  k1*  segment  having  thickness  h*  and  area  A*  (see 
Equations  (31)  and  (32))  is 


(dc) 


(87) 


Ep 

which  is  found  using  standard  analytical  methods.  The  deflection  for  the  entire  tenth  is  then 
simply  the  sum  of  the  deflection  of  each  individual  segment: 


6C  = 


JL  L 

E,  ‘ 


(88a) 
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Figure  58.  Individual  Touth  Loaded  in  Com- 
pression 


Th  e compliance  Qc  is  then  found  to  be 


Qc 


is_=_Ly  h. 

P E„  TT'  Ak 


(88b) 


Since  the  tooth  for  this  mode  of  deflection  is  in  direct  contact  with  the  material  at  the  root 
of  the  mating  tooth,  the  compliance  due  to  Hertzian  contact  deformation  must  be  included. 
Equations  (79),  (52),  (80),  (53),  and  (81)  may  be  used  for  this  purpose.  Since  the  assumed 
contracting  cylinders  are  now  concentric,  the  larger  of  the  two  radii  in  Equation  (79)  and  (52) 
must  be  negated  to  yield  a correct  equivalent  radius  term,  and  distances  d,  are  set  equal  to  the 
tooth  height  and  half  the  depth  of  its  mating  piece.  As  before,  some  average  value  will  be  assumed 
for  load  P in  Equation  (81).  The  total  radial  compliance  is  thus 


Qr*cI!*I 


(89) 


As  before,  if  the  radii  of  the  two  mating  pieces  are  the  same,  the  Hertzian  term  of 
Equation  (89)  is  dropped. 


The  stiffness  k„  is  simply  the  inverse  of  Equation  (89): 


_ 


1 

Q Radio 


(90) 
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Splint  Coupling  Stlffntaa 

Computation  of  the  stiffness  of  a spline  coupling  involves  both  lateral  and  angular  stiffness. 
The  lateral  stiffness  term  refers  to  stiffness  along  a radial  direction  passing  through  the  shaft 
centerline,  while  the  angultr  stiffness  term  refers  to  rotational  stiffness  about  an  axiB 
perpendicular  to  the  centerline,  not  to  torsional  stiffness.  In  addition,  the  stiffnesses  are  different 
for  each  of  the  two  types  of  couplings  considered  in  this  analysis,  namely  diameter  fit  and  side  fit 
unpiloted  splines,  as  shown  in  Figure  59. 


Bearing 


Major  Diameter  Fit 
Figure  59.  Two  Types  of  Spline  Couplings 


Both  types  cf  couplings  rely  on  close  tolerances,  either  on  tooth  width  (side  fit)  or  on  the 
major  diameter  for  radial  and  angular  locations.  Since  location  is  not  provided  for  any  other 
means  (a  pilot,  for  example),  the  lateral  and  angular  stiffnesses  may  be  found  from  a knowledge 
of  the  spline  coupling  geometry  and  the  tooth  stiffnesses.  The  latter  have  already  been  derived, 
while  for  the  former  the  basic  method  employed  will  be  to  consider  each  tooth  as  an  individual 
spring  and  subject  the  teeth  to  some  defined  displacement  as  shown  in  Figure  60. 


Force  Needed  to 
Impose  Disp laeement  6 


Undef  lec  ted 
Posit  ion 


Imposed  Displacement 


Individual  Teeth 


Figure  60.  Model  of  Spline  to  Determine  Overall  Stiffness 


The  force  exerted  by  the  individual  teeth  in  the  direction  of  the  load  may  then  be  calculated 
and,  hence,  the  overall  stiffness  found.  Using  this  method,  the  variation  in  tooth  stiffness  with 
applied  load  can  be  accommodated,  either  by  computing  the  force  separately  for  each  tooth  by 
iteration  or  using  an  approximate  applied  force. 

Radial  Stlftnaaa  — Major  Dlamator  Fit 

For  teeth  loaded  on  their  outer  edge,  the  displacement  of  the  teeth  radially  inward  will  be 
found  by  using  the  model  shown  in  Figure  61 . 
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Figure  61.  Model  to  Determine  Tooth  Dis- 
placement 

Consider  the  outer  edges  of  the  external  (or  male)  spline  teeth  as  describing  a circle  of 
radius  r.  This  is  circle  number  1 of  Figure  63 . Circle  number  2 coincides  with  the  major  diameter 
of  the  internal  (femaie)  spline,  that  is,  the  diameter  at  the  roots  of  the  spline  teeth.  The  two 
circles  are  initially  coincident,  so  that  the  external  spline  tooth  tips  just  touch  the  internal  spline 
tooth  roots  (Figure  59).  Circle  number  2 is  next  givin  a downward  displacement  6.  Since,  except 
for  Hertzinn  deformation  of  the  contact  zone,  the  external  and  internal  spline  members  are 
assumed  rigid,  each  of  the  teeth  will  be  compressed  inward  by  some  amount,  depending  on  its 
location. 

The  distance  Hr  represents  the  radial  compression  inward  of  a tooth  at  ingle  0 on  the 
periphery  of  circle  number  1 . For  small  values  of  displacement,  6r  can  be  closely  approximated 
by  the  equation 

6,  a 5 sin  0 0 £0^  » (91) 

The  normal  force  generated  by  the  tooth  is  the  product  of  tooth  stiffness  k*  and  the  displacement 
it: 


P*  - k.  4,  (92) 

where  stiffness  k*  is  found  from  Equation  (90) . The  contribution  to  the  force  parallel  to  the  * axis 
(i.e.,  in  line  with  the  original  displacement)  is  then 

Pj  - PR  sin  0 (98) 
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It  order  to  find  the  total  force  along  the  tj  axis  we  must  sum  all  the  individual  tooth  forces  P„.  If 
there  are  n teeth  in  the  region  0 <0<  r,  then  the  total  force  P resulting  from  an  imposed 
deflection  4 is 


P = 


n 


IP,, 


(94) 


which,  when  Equations  (91),  (92),  and  (93)  are  substituted  into  it,  becomes 


n 

P » 4 X-  kR,  sin*  $t  0 

i-» 

The  resultant  lateral  stirtneae  kL|M  for  u major  diameter  fit  spline  is  then 

n 

kL,M  = 2Z  kRI  sin*  0, 

i-t 


(95) 


(96) 


Since  the  individual  tooth  stiffness  kR|,  can  be  a weak  function  of  applied  load,  an 
approximation  of  this  applied  load  may  be  found  by  computing  the  radial  load  on  each  tooth 
using  a modification  of  Equation  (85): 


(97) 


P,  may  then  be  inserted  into  Equation  (80)  as  an  average  force  per  unit  axial  length  (hence  the 
appearance  oft  in  the  denominator  of  Equation  (97)). 


To  reiterate,  the  normal  radial  force  on  each  tooth  is  first  computed  using  only  the  term  for 
axial  compression  of  the  body  of  the  tooth.  These  normal  forces  are  then  inserted  into  the 
Hertzian  terms  for  the  stiffness  of  the  individual  teeth.  The  lateral  stiffness  iu  then  calculated 
using  these  modified  stiffness  values. 

Latmral  Stfttnoaa  — Sid*  Fit 

The  approach  for  computing  the  lateral  stiffness  will  be  basically  similar  to  that  used  for  the 
major  diameter  fit  type.  A downward  deflection  will  be  assumed  for  the  interns)  or  female  spline 
and  the  force  contribution  of  each  tooth  in  the  direction  of  the  imposed  displacement  will  be 
found  and  then  summed  over  all  teeth.  Figure  62  shows  a typical  pair  of  teeth  in  contact.  This  will 
be  defined  as  the  i1*1  contact  point. 


The  contact  point  is  at  the  pitch  radius  and  determines  angle  8.  The  angle  4 is  the  contact 
angle  and  is  determined  by  tooth  geometry.  With  the  geometry  established,  the  outer  member  is 
given  a small  downward  deflection  i which  results  in  a normal  force  P„  generated  at  the  interface 
as  the  result  of  a total  displacement  4a)  in  the  direction  of  the  normal  force.  As  before, 
displacement  4 is  small  enough  such  that  the  original  geometry  is  not  significantly  distorted  The 
angle  ♦ normal  to  the  plane  of  contact  may  be  defined  as 


♦ - 8 + 4 


(98) 
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wiere  4 may  fc«  positive  or  negative  depending  on  whether,  for  this  geometry,  the  face  is  torque- 
loaded  or  not.  If  the  leading  or  torque  load  carrying  face  of  the  male  tooth  is  also  loaded  by  PB|, 
then  0 is  negative.  The  actual  deflection  of  the  mating  teeth  normal  to  their  contact  faces  is 

= 5 1 cos*  j (99) 


Figure  62.  Model  of  Split & Teeth  in  Contact 


The  stiffness  of  the  mating  teath,  kT,  was  given  in  Equation  (86).  It  will  be  relabeled  kTi 
because  it  may  be  different  for  each  tooth  due  to  the  Hertzian  term  nonlinearity.  The  normal 
force  FP,  is  then 

Pnl  - kT,  (100) 

The  resulting  force  along  the  line  of  6 is  then  P(  for  each  tooth,  which  can  be  seen  from 
Figure  62  to  be 


P,  = Pm  | coe  ♦(  | 


(101) 


From  Equations  (99)  and  (100)  it  is  evident  that 


Pnl  = kTl  6 1 coe  *,  | 


(102) 


and  from  Equations  (101)  and  (1C2)  we  see  that 
P,  = kTI  6 cos' 


(103) 
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The  average  normal  force  is  found,  as  described  before,  for  each  tooth  by  the  equation 


= kxi  6 1 coe  I 


AVC.i  ~ "-TI 


where,  referring  to  Equation  (62), 


T Qb,  + Qs,  + Q*i  + Qb,  + Qs,  + uuo' 

In  addition,  since  the  leading  flanks  of  the  male  teeth  carry  a torque  load,  PAvg.i  must  be  modified 
to  reflect  this  factor: 

PAvc,i  = k«  « | cos  | ± (106) 

where  T is  the  transmitted  torque,  rp  is  the  pitch  radius,  and  2n  is  the  total  number  of  spline 
teeth.  The  sign  of  the  last  term  of  Equation  (106)  is  defined  as  follows: 

+ for  — */2  < 6 < t/2 

- for  t/2  < 8 < — ir/2 

The  total  lateral  stiffness  kt,*  for  the  side  fit  coupling  is  then  found  by  summing  Equation 
(103)  over  n teeth  resisting  the  deflection  and  dividing  through  by  that  deflection. 


kL,a  = kTI  cos'  ♦ 


Angular  Stlffnaaa 

The  angular  stiffness  for  a spline  coupling  will  be  a function  of  lateral  stiffness  and  spline 
geometry.  Figure  63  shows  the  right  half  of  such  a coupling,  where  the  y axis  runs  down  the 
coupling  vertical  centerline. 

The  bodies  of  the  internal  and  external  spline  halves  are  assumed  to  be  rigid  so  that  the 
deformation  of  the  teeth  alone  supplies  the  restoring  moment.  As  shown  in  Figure  63,  we  shall 
look  at  one  half  of  a coupling  rotated  about  its  left  end  plane  (the  center  of  the  entire  coupling) 
by  a small  angle  «.  Further,  a small  amount  of  radial  clearance  c is  permitted.  Since  the  bodies 
of  both  coupling  halves  are  assumed  to  be  rigid,  for  the  male  half  to  rotate  through  angle  t it  must 
deform  the  tooth  material,  starting  at  a distance  c from  the  undeflected  position.  If  a small 
section  in  the  deformed  portion  is  studied,  it  would  appear  as  in  Figure  64. 


87 


Figure  63.  Geometry  for  Calculation  of  Angular  Stiffness 


Undeformed  Surface 
Deformed  Surface 

Figure  64.  Enlargement  of  Small  Section  of  Contact  Zone 


The  stiffness  contribution  of  a small  segment  of  length  dz  is 

k..,  =>  -jgg-  dz  (108) 

where  kL  is  the  overall  lateral  stiffness  of  the  coupling  and  2b  is  the  length  of  coupling 
engagement.  The  deflection  M.Z)  is  found  from  Figure  63  to  be 

«(Z)  - t ( z + a«)-c  (109) 

and  the  force  generated  in  the  differential  segment  is  the  product  of  Equations  (108)  and  (109): 

F(Z)  " W [ * (I  + a<)_c  ] dl  (110) 

The  moment  resulting  ffom  this  force  is 

M(Z)  - ZF(Z)  - [ a (*  + a«)-c  ] dz  (111) 
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The  moment  is  next  multiplied  by  2 and  integrated  over  the  length  of  contact. 


M(Z)  = z [ « (z  + at)-c  ] dz 

where  the  limits  of  integration  d and  e are  the  boundaries  of  the  contact  area: 

j c 
d = — -a« 
e 

e = b - *(a  + be) 

When  the  integration  is  performed  and  the  result  evaluated  at  the  limits,  we  get 

M = ir  Hr (e>  ~ d,)  + (e*  - d‘> ! 

The  angular  stiffness  kA  is  then 

M kL  f (e*  - d')  , (aF  - c)  (e*  - d1)  \ 

k^  — = v*— r~  + * » 


(112) 

(113) 

(114) 

(115) 

(116) 


Because  the  kL  term  may  be  nonlinear  due  to  Hertzian  deflection,  an  average  deflection  may  need 
to  be  assumed  in  order  to  obtain  an  average  force.  (Refer  to  the  previous  stiffness  analysis  for 
major  diameter  and  side  fit  couplings.)  This  will  be  handled  by  finding  an  “average”  deflection 
Aavo  for  use  in  calculating  the  average  force  used  in  the  Hertzian  term.  The  equation  for  thiB 
deflection  is  Equation  (109)  evaluated  at  the  average  of  the  two  end-points,  Equations  (113)  and 
(114). 


«avq  = y (b«  - »•  - C)  (117) 

For  the  case  of  major  diameter  fit  couplings,  lateral  stiffness  kL  is  given  by  Equation  (96). 
Clearance  c is  assumed  to  be  zero  and  the  resulting  calculation  is  straightforward.  For  aide  fit 
spline  couplings,  there  is  a finite  clearance  between  the  major  diameters  of  the  internal  and 
external  splines,  which  is  the  clearance  c in  Figure  63.  If  c is  small  enough  (be  <c),  the  contact 
is  only  on  the  tooth  flanks  and  Equation  (106)  may  be  used  to  describe  the  lateral  stiffness,  with 
c = 0.  If,  however,  the  angle  is  large  enough  to  cause  contact  at  the  major  diameter,  then 
Equation  (116)  must  be  applied  a second  time,  using  Equation  (96)  for  radial  stiffness  and  with 
c being  set  equal  to  the  major  diameter  clearance. 

to.  Spun*  Coupling  Damping  CoafUdanta 

Damping  coefficients  for  a spline  coupling  are  found  by  assuming  the  damping  to  consist 
primarily  of  Coulomb  damping.  The  energy  dissipated  per  cycle  by  this  mechanum  will  be 
computed  and  used  to  find  an  equivalent  viscous  damping  coefficient  for  use  in  the  stability 
calculation.  Figure  66  shows  an  elementary  cpring-maas-damper  system  where  the  damping  is 
provided  by  a rigidly  connected  Coulomb  damper. 


Rh.  1 


y/)/ /'//}// 

Figure  66.  Elementary 
Spring  — Maas  - 
Coulomb  Damper 
System 

The  differential  equation  cf  motion  for  this  system  is 
mi  + ki  i P(  = 0 


(118) 


where  the  frictional  force  Ff  is  constant  in  magnitude  and  always  in  a direction  opposed  to  the 
velocity.  An  approximate  solution  to  this  equation  is  based  on  the  equivalence  of  energy 
dissipation  between  this  and  a viscous  damped  system.  The  energy  Ec  dissipated  per  cycle  for  a 
Coulomb  damped  system  is  simply 


Ec  - 4Fr  «„ 


(119) 


where  5 is  the  maximum  zero-to-peak  amplitude  of  the  system.  For  a viscous  damped  system,  the 
energy  dissipated  per  cycle  Ev  is 


E,  = xCw6J 


(120) 


where  C is  the  damping  coefficient,  w is  the  angular  frequency,  and  S„  is  as  previously  defined.  If 
we  wish  to  define  a viscous  damped  system  equivalent  tc  the  >ne  in  Figure  65,  Equations  (119) 
and  1 120)  are  equated  and  solved  for  C,„,  the  equivalent  viscous  damning  coefficient: 


C- 


xwbo 


(121) 


A coefficient  has  thus  been  found  which,  for  a specified  oscillation  amplitude  and  frequency, 
causes  the  same  energy  to  be  withdrawn  using  a viscous  mechanism  as  the  original  system  did 
using  sliding  friction.  The  ,wo  systems  would  then  be  expected  to  exhibit  similar  behavior  under 
these  conditions. 

The  spline  coupling  damping  coefficients  will  be  treated  in  the  above  manner,  where  will 
be  the  relative  lateral  or  angular  displacement  of  the  coupling  pieces  and  F(  the  total  frictional 
force  or  moment  resulting  from  b0.  Since  instability  is  expected  to  occur  at  the  first  natural 
frequency,  w is  the  difference  between  the  running  speed  (or  speed  of  mtetest)  and  the  first 
natural  frequency  where  the  former  is  assumed  always  greater  than  the  latter. 


The  friction  force  may  be  computed  using  the  equation 
Ft  » f Fn 


(122) 
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where  f is  the  coefficient  of  friction  and  F„  the  normal  force.  The  overall  coupling  friction  force 
will  he  found  by  computing  the  normal  force  on  each  pair  of  sliding  surfaces,  applying  an  assumed 
or  computed  coefficient  of  friction,  and  summing  the  friction  force  components  for  all  the 
interfaces. 


If  the  coupling  teeth  ere  evenly  spaced,  the  force  in  the  circumferential  direction  FT  due  to 
applied  torque  T is 


Ft  = 


T 

r„n 


(123) 


where  n is  the  number  of  teeth  and  rp  is  the  pitch  radius.  Figure  66  illustrates  the  geometry  for 
computation  of  the  normal  force  FeTfor  a single  tooth.  From  the  figure,  it  can  be  seen  that 


F 


n.T  ~ 


_2i_ 

COQ0 


(124) 


Hence,  the  normal  force  on  the  tooth  flank  due  to  torque  is  found  from  Equations  (123)  and  (124) 


to  be 

F,t  = 


r^\ 

1 


rpn  coa  0 


(126) 


The  friction  force  for  that  toot.h  may  be  found  from  Equations  (122)  and  (126)  to  be 


F 


n.T  = 


fT 

r„n  cos  <t> 


(126) 


The  energy  dissipated  per  cycle  for  the  spline  coupling  will  now  be  computed.  Since  the 
oscillation  caused  by  an  instability  will  ..ot  be  laterally  planar  but  will  rather  be  a whirl 
phenomenon,  all  contacting  interfaces  are  fully  cycled  for  each  cycle  of  spline  motion.  The  energy 
dissipated  per  cycle  is  then  the  energy  dissipated  at  one  interfacer-for  a completed  cycle,  multipled 
by  the  number  of  such  interfaces.  This  method,  of  course,  assumes  circular  nonsynchronous 
whirl. 
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If  a radial  displacement  5„  estate,  then  the  interface  moves  a distance  6t  such  that 


Equation  (119)  describes  the  energy  dissipated  per  cycle  for  a general  friction  damped 
system.  When  Equations  (126)  and  (127)  are  substituted  into  (119),  the  result  is 


which  is  the  energy  dissipated  by  one  interface  over  a full  cycle.  For  n interfaces,  the  energy 
dissipated  Ec|R  in  radial  motion  per  cycle  is 


rpcos‘  0 


Note  that  there  is  no  dependence  on  the  number  of  spline  teeth  in  contact  (using  this 
particular  model). 

For  angular  motion,  the  energy  for  V*  cycle  is  found  by  assuming  an  angular  displacement 
of  a,„  as  shown  in  Figure  67. 


Spline- 

£ 


Dif  ferent ia  1 
Contact  Area 


Figure  67.  Geometry  for  Frictional  Moment 


The  friction  force  dFf  for  a small  section  of  the  contact  area  is 


dF,  « dz  (130) 

where  2b  is  the  length  of  the  spline  and  dz  is  the  width  of  the  contact  area. 

The  distance  through  which  this  force  moves  is  za«  /cos  *.  The  energy  dissipated  by  one 
interface  is  then 


E - 2 / d* 
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where  the  factor  of  2 accounts  for  the  ether  half  of  the  interface  (negative  z).  The  energy 
dissipated  per  cycle  is  the  above  integral,  evaluated  at  its  limits,  and  multiplied  by  a factor  of 
four: 


E„i  = 


2F  rctob 
coe  <p 


(132) 


When  Equation  (126)  is  substituted  into  (132)  and  the  result  multipled  by  the  number  of 
interfaces,  the  energy  dissipated  by  the  spline  per  cycle  of  angular  motion  (conical)  results: 


E 


c.A 


2f  Tbop 
rp  coe 4 


(133) 


Equivalent  viscous  damping  coefficients,  b and  0,  for  lateral  and  angular  motion,  respectively, 
can  be  found  using  Equation  (121).  Since  b and  (t  cue  for  one  axis  only,  Equations  (129)  and  (133), 
which  cover  complete  motion,  must  be  divided  by  a factor  of  two  before  being  substituted  into 
Equation  (121).  The  net  results  of  these  substitutions  are 


b 


2fT 

TpXwio  CO6X0 


0 


fTb 

rpx-ua0  cos*$ 


(134) 

(136) 


The  two  values  above  may  be  used  directly  as  the  equivalent  viscous  coefficients  which  simulate 
for  given  conditions  the  friction  damping  present  in  the  coupling. 

A description  of  the  complete  computer  program  (Report  No.  FR-10601)  was  provided  to  the 
Applied  Technology  Laboratory  in  October  1978.  A User's  Manual  for  this  program  appears  in 
Appendix  B. 


c.  Strategy  for  Stability  Calculation a 

The  following  is  a suggestion  for  the  strategy  to  be  used  in  performing  a stability  calculation 
for  a particular  rotor  bearing-coupling  system: 

1.  Model  the  rotor-bearing  system  using  standard  procedures.  Define  the  spline 
coupling  geometry  and  torque. 

2.  Assume  a running  speed  and  vertical  and  lateral  whirl  amplitudes.  Since 
amplitude  appears  in  the  equivalent  viscous  damper  expression  as  a 
denominator,  the  damping  becomes  undefined  for  zero  whirl  amplitude.  The 
following  two  procedures  for  estimating  the  amplitude  are  suggested: 

a.  If  the  misalignment  is  small,  the  coupling  may  lock  up  for 
small  whirl  amplitudes,  accommodating  small  displace- 
ments through  tooth  deflection.  The  maximum  friction 
force  at  a tooth  interface  is  given  by  Equation  (126),  where 
f is  now  the  coefficient  of  static  friction.  This  represents  the 
largest  force  transmittabie  by  a pair  of  teeth  in  contact  with 
only  elastic  tooth  deflection.  The  deflection  of  the  mating 
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teeth  along  the  line  of  force  is  given  by  the  sum  of 
Equations  (68),  (64)  and  (71)  for  each  tooth.  If  the 
amplitude  is  assumed  to  be  slightly  larger,  it  cannot  be 
accommodated  by  elastic  deflection,  sliding  occurs  with  the 
accompanying  internal  friction,  and  instability  is  a possible 
result  if  external  damping  is  not  large  enough.  We  have 
thus  established  a “floor”  for  our  assumed  amplitude. 


b.  When  misalignment  is  large,  sliding  must  already  occur  in 
order  for  the  coupling  to  accommodate  it.  The  coupling, 
therefore,  cannot  lock  up,  and  it  is  possible  for  any  whirl 
amplitude  to  generate  internal  sliding  friction.  For  this 
case,  it  appears  best  to  assume  a whirl  amplitude  for  the 
shaft  system  in  the  first  mode,  as  it  is  this  mode  which  is 
expected  to  be  excited.  The  coupling  lateral  and  angular 
amplitude  may  be  easily  formed  from  the  assumed  de- 
flection and  the  mode  shape. 


3.  Using  the  methods  outlined  in  Sections  V.A.2.a  and  V.A.2.b,  calculate 
coupling  stiffnesses  k and  H and  damping  coefficients  b and  /J.  The 
coefficient  of  friction  must  be  estimated  in  order  to  find  the  latter  two 
coefficients.  It  is  assumed  that  the  coupling  is  lubricated  and,  furthermore, 
enough  lubricant  gets  between  the  contacting  faces  to  prevent  metal-to- 
metal  contact.  The  relatively  low  sliding  velocities  (a  few  tens  of  centimeters 
per  second),  the  high  loads,  and  the  oscillatory  motion  suggest  that  boundary 
lubrication  (as  opposed  to  hydrodynamic  or  elastohydrodynamic)  occurs. 
This  suggests  that  friction  coefficients  will  generally  be  not  too  far  below 
those  seen  in  metal-to-metal  contact  and  certainly  greater  than  those 
resulting  from  a film  generated  at  higher  velocities  or  through  a 
rolling/sliding  mechanism.  Reference  17  contains  a plot  of  coefficient  of 
friction  against  sliding  speed  for“surfaces  capable  of  hydrodynamic  lubri- 
cation." In  the  boundary  lubrication  region,  the  coefficient  of  friction  lies  in 
the  approximate  range  from  0.08  to  0.10.  Reference  18  gives  an  empirical 
equation  for  coeificient  of  friction  for  the  rolling-sliding  case,  derived  from 
roller  test  data: 


f 


0.0127  log,. 


3.17  x 10* 

(fa^JL  vT* 

w 1 


(136) 


where  m°  is  the  lubricant  viscosity  (cp),  w is  the  load  (tb/in.  of  face  width),  and 
V,  and  VT  are  sliding  and  sum  velocities  (in./sec)  of  the  two  surfaces, 
respectively.  Realising  that  we  have  pure  sliding  and  assuming  V,  VT*  can  be 


17  Ribtnowits,  E. “Frictional  and  Wear  Malarial*, "John  Wiley  and  Sent,  New  Yort:,  1986. 

18  Benedict,  G.  H„  and  Kelly,  B.  W “Instantaneous  Coefficient*  of  Gear  Tooth  Friction,' "ASUS  Transaction*  4,  69-79, 
1961 
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replaced  with  the  cube  of  the  sliding  velocity  V,*,  we  can  obtain  an  equation 
for  further  guidance  on  the  value  of  the  friction  coefficient: 


f = 0.0127  logm 


3.17  X 10*  ” 

JUsl  V • 

w * _ 


(137) 


which  is  valid  for  values  of  the  friction  coefficient  from  about  0.005  to  0.08. 

Boylan  (Reference  19)  performed  tests  on  gear  couplings  and  measured 
coefficients  of  friction  as  a result  of  misalignment  (replotted  in  Reference  20 
vs  sliding  velocity).  Coefficients  of  friction  ranged  from  0.15  to  0.25  for  low 
velocity  to  a low  of  0.05  to  0.10  at  about  60  cm/sec  sliding  velocity  and  rose 
above  that  velocity.  These  tests  were  done  on  a test  rig  capable  of  rotational 
speeds  to  9000  rpm,  and  the  friction  was  apparently  very  difficult  to  pin 
down. 

The  velocity  may  be  calculated  for  lateral  motion  by  the  equation 

Vmax.l  = 60  (138) 

and  for  angular  motion  the  “average”  maximum  velocity,  that  is  the  maximum 
velocity  at  a point  halfway  between  the  center  and  edge,  is 

Vmax.a  a (139) 

An  average  velocity  for  both  cases  would  then  be  the  maximum  divided  by  \f~2-. 

vavo,l  = a (j4°) 

v*VQ.A  = 77T  n (141) 

Since  angular  and  lateral  misalignment  probably  will  occur  along  the  same 
plane,  velocity  totals  would  be  the  sums  of  Equations  (138)  and  ( 139)  or  ( 140) 
and  (141).  Equations  (138)  and  (139)  would  be  used  for  the  friction- 
correlation  of  References  19  and  20,  while  Equations  (140)  and  (141)  appear 
to  apply  for  the  cases  of  References  17  and  18. 

4.  Having  all  necessary  coefficients  in  hand,  the  damped  eigenvalues  of  the 
system  can  now  be  computed.  This  should  be  done  for  several  speeds,  loads, 
and  misalignments  in  the  operating  range  to  obtain  a clear  stability  map 
(plot  of  the  eigenfrequencies  versus  the  varied  parameters)  of  the  system 
being  studied. 


19  Boyland,  W. .“Marin*  Application  of  Dental  Couplings,"  Soc.  Nav.  Arch,  and  Nav.  Engr*.  Marin*  Power  Plant 
Symposium.  Paper  2S,  May  1966. 

20  Cr**M,  A.  B. "Design  Principim  and  Lubrication  of  Gear  Coupling* . 'international  Conference  on  Flexible  Coupling* 
for  High  Power*  and  Speeds,  University  Df  Suaaei.  England.  June  2ft  to  July  1.  1977,  Paper  Bl. 
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8.  EXPERIMENTAL  RIG  DESIGN 


The  nonsynchronous  rig  configuration  is  shown  in  Figures  68  and  69.  It  consists  of  a damped 
thrust  bearing  simulator  rig,  which  induces  a turbine  shaft  supported  on  two  bearings.  The  rear 
end  of  the  shaft  is  supported  on  a duplex  pair  ball  bearing  and  a squeeze  film  damper  through  a 
flexible  bearing  support.  The  front  of  the  shaft  has  been  modified  to  accommodate  the 
interchangeable  test  splines  which  are  supported  by  a separate  set  of  duplex  ball  bearings.  The 
rig  was  driven  by  a turbine  which  is  attached  at  the  rear  of  the  existing  rotor  dynamics  rig  shaft. 
High  pressure  steam  flowed  through  the  turbine  to  develop  the  torque  necessary  to  produce 
nonsynchronous  whirl.  The  torque  was  transmitted  through  the  spline  and  absorbed  via  a water 
brake.  The  transmitted  torque  was  measured  at  the  water  brake  with  a torque  arm  mounted  on 
a force  scale. 

The  turbine  end  ball  bearings  are  axially  preloaded  against  a second  pair  of  bearings 
through  a connecting  rod  to  prevent  bearing  skid.  The  test  spline  stub  shaft  is  also  supported  on 
preloaded  duplex  bearings  as  shown  in  Figure  70.  To  allow  for  independent  control  of  the  spline 
lubrication,  a separate  supply  line  is  provided  to  the  spline  area  (item  1 in  Figure  70).  In  addition, 
a seal  cup  is  provided  so  that  the  bearing  oil  will  not  lubricate  the  spline  teeth  during  the  “dry  ” 
spline  portion  of  the  test  (item  2 in  Figure  70).  The  entire  rig  lubrication  schematic  is  shown  in 
Figure  71. 

Figures  72  and  73  are  photographs  of  the  rig  as  installed  in  the  test  stand.  The  former  shows 
an  overall  view  of  the  turbine  end  of  the  rig  and  the  thrust  bearing  housing.  The  latter  shows  an 
overall  view  of  the  water  brake  end  of  the  rig  and  the  test  spline  housing. 


C.  RIG  TEST 
1.  Test  Plan 

The  objective  of  the  nonsynchronous  response  rig  test  was  to  generate  experimental  data  to 
verify  the  analytical  predictions  of  the  nonsynchronous  whirl  phenomenon  induced  by  spline 
friction.  Parameters  such  as  damper  oil  supply,  type  of  spline,  spline  clearance  fit,  amount  and 
type  of  spline  lubricant,  and  their  effects  on  spline  induced  vibration  was  studied.  The  test  data 
was  correlated  with  the  prediction  of  the  spline  friction  analytical  model  developed  in  Section 
V.A. 


The  instrumentation  is  depicted  in  Figure  74.  Shaft  response  was  monitored  by  means  cf 
vertical  and  horizontal  accelerometers  located  at  the  turbine  end  bearing  support.  In  addition, 
two  proximity  probes  at  the  midspan  of  the  shaft  measured  midspan  shaft  deflection.  The  signals 
from  the  vertical  accelerometer  and  proximity  probe  were  processed  through  a real-time  analyzer 
for  on-line  data  evaluation.  This  provided  continuous  monitoring  of  all  nonsynchronous  rotor 
vibrations.  In  addition  to  the  real-time  analyzer,  the  vertical  and  horizontal  proximity  signals 
were  fed  through  a tracking  filter.  The  filtered  and  unfiltered  signals  were  plotted  to  isolate  the 
nonsynchronous  vibration  components.  All  signals  were  recorded  on  magnetic  tape.  A photograph 
of  the  test  stand  instrumentation  setup  is  presented  in  Figure  75. 

The  test  program  is  outlined  in  Figure  76.  The  test  hardware  consisted  of  four  different 
splines:  side  fit  and  major  diameter  fit  with  two  different  clearances  each  (designed  according  to 
ANSI  spec  B92.1  - 1970). 
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Figure  69.  Detail  of  Drive  Turbine  Assembly 


The  test  program  called  for  sis:  baseline  runs.  Each  run  consisted  of  a slow  acceleration 
through  the  rig  operating  speed  range  with  a “dry"  spline  and  damper  as  a baseline.  Once  a 
threshold  of  instability  was  establishsd,  spline  lubricant  was  allowed  to  flow  at  three  different 
rates.  The  effect  of  each  flow  rate  on  the  instability  speed  was  recorded.  Once  the  threshold  had 
been  reached  at  each  spline  flow  rate,  damper  lubricnnt  was  allowed  to  flow  to  the  squeeze  Him 
damper.  The  damper  was  utilized  to  suppress  spline-induced  instabilities  throughout  the  rig 
speed  range. 


i 
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Lubrication  System 


Figure  73.  Nonsynchronous  Rotor  Dynamics  Rig  D-8  Stand, Overall  View,  Front 


Freq 

Figure  74.  Nonsynchronous  Rig  Instrumentation  Schematic  Diagram 


figure  75.  Nonsyriclironous  Rotor  Dynamics  Rig  D-8  Stand  Instrumentation  Setup 


Spline  ANSI 
B92-1 

Side  Fit 

Major  Die 

imeter  Fit 

1970 

CL  1* 

CL2 

CL  1* 

CL  2 

Run  Mode 

0.0036 

0.0065, 

0.0036 

0.0012. 

0.00 15L 

0.0035~ 

0.0026 

0 0002 

Dry  Spline  Baseline 

X 

X 

X 

X 

Spline  0.15Appm 

X 

X 

X 

X 

Lubricant  020Bppm 

X 

X 

X 

X 

Q 

Flow  rate  n _ 

0.25  ppm 

X 

X 

X 

X 

Active  Damper 

X 

X 

X 

X 

•These  Columns  Were  Repeated  for  Type  III  Lubricant 
Mar.  Rig  Speed  7000  rpm 
Max  Rig  Torque  4500  in.-ib 

Figure  76.  Nonsync hronous  Test  Program 
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2.  T«*t  Result* 


The  nonsynchronous  whirl  test  rig  results  are  summarized  in  Figure  77.  Shown  in  the  same 
form  as  the  test  program  (Figure  76),  this  Figure  illustrates  which  of  the  test  rig  configurations 
experienced  spline-friction-induced  rotor  instabilities.  The  side  fit  spline  was  unstable  for  both 
the  loose  and  tight  tooth  clearance  configuration  when  unlubricated.  In  all  test  runs  (both  side  fit 
and  major  diameter  fit  splines),  when  spline  lubrication  was  allowed  the  rotor  instability  was 
suppressed.  The  major  diameter  fit  splines  were  always  stable;  even  when  unlubricated  the  spline 
lubricant  flow  rate  had  no  effect  on  stability.  To  determine  the  effect  of  squeeze  dim  damping  on 
rotor  instability,  the  tight  clearance  side  fit  spline  was  run  unlubricated  with  the  rig  squeeze  film 
damper  activated.  The  result  was  that  the  damper  also  acted  to  suppress  the  unstable 
nonsynchronous  whirl. 


Spline 

Side  Fit 

Maior  Diameter  Fit 

Configuration 

Run  Mode 

Loose 

Tight 

Loose 

Tight 

“Dry” 

Spline 

U 

U 

S 

S 

"Wet" 

Spline 

S 

S 

S 

S 

“Dry"  Spline 
“Activated"  Damper 

S 

U ~ Unstable  Rotor  Vibrations  (Nonsynchronous) 
S = Stable  Rotor  Vibrations  (Synchronous) 


Figure  77.  Summary  of  Nonsynchronous  Rig  Test  Results 


The  unfiltered  midshaft  vibration  levels  for  the  tight  clearance  side  fit  spline  are  presented 
in  Figure  78.  The  response  plot  shows  e rotor  bounce  resonance  at  about  2400  rpm,  and  at  about 
4200  rpm  the  overall  vibration  level  climbs  very  rapidly  for  the  unlubricated  spline,  but  there  is 
no  increase  for  the  lubricated  spline.  Rotor  instability  occurred  at  that  poii.t,  and  this  is  verified 
by  the  frequency  response  plot  of  Figure  79.  Here,  the  frequency  components  of  the  4200  rpm  are 
shown.  Note  that  the  largest  component  of  response  is  at  64%  of  the  shaft  speed  (0.64E).  Also,  a 
small  2E  vibration,  which  can  be  associated  with  a slight  coupling  misalignment,  can  be 
observed. 

For  comparison,  Figure  80  shows  the  frequency  response  plot  of  lubricated  side  fit  spline  test 
at  7000  rpm.  Note  that  no  subsynchronous  frequencies  are  present  and  only  a small  2E 
component  can  be  observed.  All  of  the  runs  summarized  in  Figure  77  had  the  same  response 
characteristics  as  those  shown. 
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Figure  78.  Tight  Clearance  Side  Fit  Spline  Rotor  Response  (Lubricated  and 
Unlubricaied) 
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Frequency  Hz 

Figure  80.  Response  Amplitude  vs  Frequency  at  7(00  rprr  \ pub  icated  Side  Fit 
Spline) 


In  addition  to  the  planned  tests,  one  run  • as  made  to  verify  the  effectiveness  of  the  rig’s 
squeeze  film  damper  to  reduce  the  steady -state  unbalance  response  amplitude.  With  very  low 
water  brake  load,  t^e  ~ig  shaft  speed  was  increased  from  0 to  7000  rpm.  Figure  81  shows  the 
midshaft  vibratir  n response  with  and  without  the  squeeze  film  damper  activated.  Note  that  at 
least  a 50'*.’  reduction  in  rotor  response  occurred  when  the  rig  was  damped.  The  increase  in 
damper  oil  supply  pressure  had  very  little  effect  on  rig  response. 


D.  DATA  CORHELaTION  WITH  ANALYTICAL  MODEL 

The  analytics'  model  used  to  predict  the  debility  characteristics  of  the  nonsynchronous 
whirl  rig  was  the  spline  coupling  friction  model  developed  in  Task  VIII  of  this  program  and  used 
with  a rotor  syt  tern  stability  m<xlel  similar  to  that  of  Reference  13.  The  rotor  and  spline  physical 
dimensions  were  input  to  the  stability  model,  and  the  output  was  complex  eigenvalues  and 
eigenvectorc.  The  imaginary  part  of  the  eigenvalue  gives  the  frequency  of  vibration  and  the  real 
part  gives  the  exponential  growth  factor.  lr  positive,  the  vioration  is  unstable;  if  negative,  it  is 
stable.  The  eigenvectors  give  the  mode  shape  of  the  whirling  shaft.  A parametric  study  was 
performed  to  find  the  effects  of  changing  various  rig  and  coupling  parameters  on  the  rig  stability 
characteristics. 

The  rig  was  modeled  for  use  in  the  rot'  r dynamics  programs  with  a 12-station,  11 -section 
breakup.  The  model  was  run  to  determine  its  undamped  critical  speeds  with  the  assumption  that 
the  spline  coupling  is  very  stiff.  The  fust,  four  undamped  critical  speeds  were  found  to  occur  at 
2376,  5092,  9354,  and  14,885  rpm.  The  mode  shapes  for  the  first  two  speeds  are  shown  in 
Figures  82  and  83. 
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Critical  Speed  ■ 5092  rpm 
Figure  83.  Second  Critical  Speed  Nonsynchronous  Rig 


For  subsequent  stability  runs,  the  model  was  modified  by  inserting  the  spline  coupling 
stiffness  and  damping  coefficients  at  the  spline  end  of  the  shaft  and  inserting  the  damper  bearing 
coefficient  of  damping  at  the  squeeze  film  damper  location.  Basically,  the  model  describes  the 
spline  coupling  in  terms  of  four  coefficients: 

1.  Radial  stiffness,  k (lb/in.) 

2.  Angular  stiffness,  « (lb/in. radian) 

3.  Radial  damping,  b (lb-sec/in.) 

4.  Angular  damping,  0 (in.-lb-sec/radian) 

These  values  are  then  used  with  the  rotor  model  as  input  to  the  rotor-bearing  system  stability 
program. 

The  angular  stiffness  of  the  coupling  varies  with  the  angular  whirl  amplit  ’ he  nplitude 
of  the  external  spline  coupling  half  within  the  internal  half).  The  two  dampin'  .onts  are 

sensitive  to  the  following  quantities: 

• Running  speed 

• Coupling  angular  whirl  amplitude  (external  spline  relative  to  internal) 

• Rotor  system  first  critical  speed  (at  which  instability  is  expected  to  occur) 

• Coupling  radial  whirl  amplitude  (external  spline  relative  to  internal) 

• Coefficient  of  friction  of  spline  coupling 

• Transmitted  torque. 
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Because  the  Coulomb  friction  damping  actually  present  in  the  coupling  has  been  modeled 
with  equivalent  viscous  damping  coefficients,  the  resulting  coefficients  b and  $ very  with 
amplitude,  relative  speed  (running  speed  minus  first  critical  speed),  pressure  angle  at  the  pitch 
diameter  (a  function  of  spline  geometry),  and  all  the  parameters  previously  listed.  A spline 
analysis  using  this  method  therefore  takes  the  form  of  a parametric  study. 

The  spline  coupling  model  was  run  with  the  following  assumptions: 

1.  The  coefficient  of  friction  is  0.15. 

2.  The  transmitted  torque  is  2000  in. -lb,  which  is  about  half  of  the  full  load 
torque  of  4052  in.- lb. 

3.  The  radial  amplitude  is  equal  to  the  radial  clearance  plus  a lateral  deflection 
of  2 X 10~*  in. 

4.  The  angular  amplitude  is  equal  to  the  angle  e defined  by  the  radial  clearance 
and  the  coupling  length  plus  a nominal  coupling  angular  deflection. 

^ . . 2C 

e £ tan~*  ~j~ 


where  C is  the  clearance  and  I is  the  coupling  length.  The  coupling  angu'ar 
deflection  was  taken  to  be  0.0076  deg.  Values  for  both  lateral  and  angular 
deflections  were  found  by  estimating  the  intercoupling  forces  and,  for 
calculated  stiffnesses,  finding  the  resulting  deflections. 

5.  The  first  critical  was  taken  from  the  undamped  analysis  to  be  2400  rpm. 

6.  The  rotor  was  running  2500  rpm  above  its  first  critical  speed,  giving  an 
operating  speed  of  4900  rpm. 

7.  The  given  fillet  radius  of  0.027  in.  does  not  define  a fillet  which  extends 
uniformly  from  the  root  of  the  tooth  to  the  form  diameter.  The  smallest  value 
for  this  radius,  assuming  a circular  fillet,  is  0.040  in.,  and  that  was  used. 

The  final  models  for  the  four  test  splines  are  given  in  Table  4.  These  served  directly  as  input 
to  the  spline  coupling  computer  program. 

The  analysis  of  the  test  rig  was  performed  by  choosing  and  analyzing  a nominal  coupling 
configuration  and  then  analyzing  perturbations  of  this  case.  For  this  study,  the  loose,  side  fit 
spline  coupling  was  defined  as  the  nominal  case,  which  is  described  in  the  third  column  in  Table 
4.  Perturbations  of  this  case  consisted  of  changes  to  the  spline  coupling  coefficients  and  changes 
to  the  test  rig  geometry  or  operating  conditions. 
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TABLE  4 

MATHEMATICAL  MODELS  OF  FOUR  TEST  RIG  SPLINE  COUPLINGS 


Loose  Major 
Diameter  Fit 

Tight  Major 
Diameter  Fit 

Loose 

Side  Fit 

Tight 

Side  Fit 

Operating  Speed  (rpm) 

4900 

4900 

4900 

4900 

Number  of  Tooth  Pairs 

36 

36 

36 

36 

Number  of  Segments,  Ext. 

10 

10 

10 

10 

Number  of  Segments,  Int. 

10 

10 

10 

10 

Pressure  Angle  (deg) 

30 

30 

30 

30 

Face  Width  (in.) 

1.75 

1.75 

1.75 

1.75 

Young’s  Modulus  (Ib/in.*) 

30  X 10* 

30  X 10* 

30  X W 

30  X 10* 

Poisson's  Ratio 

0.3 

0.3 

0.3 

0.3 

Contact  Patch  Width  (in.) 

0.1 

0.1 

0.1 

0.1 

Angle  of  Tilt  (deg) 

0.1025 

0.0239 

0.3448 

0.2446 

First  Critical  (rpm) 

2400 

2400 

2400 

2400 

Radial  Amplitude  (in.) 

1.522  X 10-* 

3.52  X 10-4 

5.252  X 10- r' 

3.752  X 10- • 

Coefficient  of  Friction 

0.15 

0.15 

0.15 

0.15 

Torque  (in. -lb) 

2000 

2000 

2000 

2000 

Radial  Clearance  (in.) 

1.55  X 10-* 

3.5  X 10-4 

5.25  X 10- * 

3.75  X 10-* 

Shear  Modulus  (Ib/in.*) 

11.58  X 10* 

11.58  X 10* 

11.58  X 10* 

11.58  X 10* 

Space  Width  (in.) 

0.1344 

0.1344 

0.1344 

0.1344 

Tooth  Thickness,  Ext.  (in.) 

0.1258 

0.1258 

0.1239 

0.1239 

Pitch  Diameter 

3.0 

3.0 

3.0 

3.0 

Minor  Diameter,  Ext.  (in.) 

2.871 

2.871 

2.83 

2.83 

Form  Diameter,  Ext.  (in.) 

2.910 

2.910 

2.91 

2.91 

Fillet  Radius,  Ext.  (in.) 

0.040 

0.040 

0.040 

0.040 

Major  Diameter,  lnt.  (in.) 

3.084 

3.084 

3.170 

3.170 

Form  Diameter,  Int.  (in.) 

3.069 

3.069 

3.090 

3.090 

Fillet  Radius,  Int.  (in.) 

0.040 

0.040 

0.040 

0.040 

In  addition  to  varying  coupling  parameters,  changes  in  the  test  rig  were  explored  for  their 
effect  on  stability.  The  parameters  which  were  varied,  including  those  of  both  coupling  and  rig, 

are: 


• Angular  stiffness  (*)  of  the  spline  coupling 

• Damping  coefficients  (b,  ff)  of  the  coupling,  varied  together 

• Stiffness  of  the  damper  bearing 

• Damping  coefficient  of  the  damper  bearing 

• Running  epeed 

• Diameter  of  the  spline  coupling  “neck”  (the  reduced  diameter  transition 
piece  from  the  main  shaft  to  the  external  spline  coupling  half) . 


The  result#  of  this  parametric  study  are  presented  in  Table  6.  The  table  shows,  for  each  run 
of  the  stability  program,  which  parameter  was  varied,  the  new  value,  and  the  resulting  damped 
natural  frequency  and  logarithmic  decrement.  Changca  in  the  coupling  stiffness  and  damping  are 
given  as  the  factor  used  for  multiplying  or  dividing  the  nominal  value;  the  others  give  the  actual 
number  used  in  place  of  tha  nominal.  These  nominal  values  t.re  given  for  reference  in  Table  6. 

Results  of  the  computer  runs  made  tc  describe  spline  coupling  coefficients  and  stability 
analysis  results  are  given  in  Tables  4 and  5.  These  results  of  the  parametric  study  an;  also 
presented  graphically  in  Figures  84  through  88.  Figure  84  shows  the  variation  of  system  log 
decrement  with  the  damping  coefficient  of  the  damper  bearing  and  illustrates  just  how  small,  for 
the  otherwise  nominal  rig  configuration,  the  damping  must  gee  to  bring  about  instability. 
Figure  86  shows  the  variation  in  log  dec  with  changes  in  the  anguiar  stiffness  of  the  spline 
coupling.  For  this  plot  the  bearing  has  been  hard-mounted  und  damping* reduced  to  a value  of 
2 !b-sec/in.  Note  that  reducing  the  angular  stiffness  has  a significant  effect  on  stability.  Stiffness 
could  be  reduced  by  shortening  the  length  of  engagement  or  crowning  (he  coupling  The  variation 
of  system  log  dec  with  spline  coupling  damping  coefficient  is  shown  in  Figure  86.  The  log  dec  is 
veiy  insensitive  to  coupling  damping  until  unrealistically  high  damping  values  are  reached. 
Figure  87  shows  that,  for  a low  value  of  bearing  damping  and  the  spline  damping  increased  by  a 
factor  of  100,  increasing  the  damper  bearing  support  stiffness  causes  the  rig  to  go  marginally 
unstable.  Finally,  Figure  88  shows  the  change  in  the  natural  frequency  with  spline  angular 
stiffness  for  a lightly  damped,  stiff  damper  bearing. 

The  predicted  first  critical  speed  of  the  nonsynchronous  whirl  rig  was  2376  rpm  and,  as 
shown  in  Figure  89,  this  calculated  value  is  in  very  close  agreement,  with  the  observed  critical 
speed.  Under  nominal  operating  conditions  the  test  rig  was  predicted  to  be  marginally  unstable 
if  very  low  damping  was  present  at  the  bearing  supports.  The  experimental  results  showed  that 
the  rotor  was  slightly  unstable  with  the  side  fit  spline  unlubricated,  but  by  reducing  the  spline 
friction  forces  with  lubricating  oil  the  instability  was  suppressed.  Also,  as  predicted  with  the  rotor 
stability  model,  squeeze  film  damping  the  rig’s  thrust  bearing  eliminated  the  nonsynchronous 
whirl.  Overall,  there  was  close  agreement  between  the  predicted  and  measured  response  of  the 
nonsynchronous  whirl  rig.  A further  description  of  the  spline  friction  stability  analysis  can  be 
found  in  References  21  and  22. 


21  Tecta,  J.  A.,  "Stability  Modal  and  Analytic  for  Spline  Coupling  Induced  Noneynchronoue  Rotor  Vibrations,” 
Mechanical  Technology  Incorporated  Report:  MTI78TR9,  October  1977. 

22  Tecta,  .T  A.  "Stability  Analyeia  of  a Spline  Coupling  Teet  Pig,”  Mechanical  Technology  Incorporated  Report: 
MTI-78TR78,  May  1973. 
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System  Log  Decrement 


TABLE  6 

NOMINAL  VALUES  FOR  PARAMETERS 


Parameter 

Spline  Angular  Stiffness 

Spline  Damping,  Radial 
Spline  Damping,  Angular 
Damper  Bearing  Stiffness 
Damper  Bearing  Damping 
Running  Speed 
Spline  Neck  Diameter 


Nominal  Value 

3.3119  X 10T  Ib-in./radian 
123.47  lb-sec/in. 

47.143  in.-lb-sec/radian 
86,000  Ib/in. 

50  lb -sec/in. 

4900  rpm 
1.78  in. 


Figure  84.  Variation  of  System  Log  Decrement  With  Bearing  Damping 
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Rotor  Speed  - rpm 

Figure  89.  Ncnsynchronous  Whirl  Rig  Critical  Speed  Observed  and  Predicted 
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SECTION  VI 


ROTOR  DESIGN  OPTIMIZATION  (TASK  X) 


At  the  conclusion  of  Tasks  VI  through  IX  the  effects  of  five  fundamental  design  variables 
(support  stiffness,  support  damping,  bumper  clearance,  spline  geometry,  and  spline  lubrication) 
on  the  dynamic  characteristics  of  damped  flexible  rotor  systems  had  been  evaluated. 
Experimentally  verified  analytical  methods  had  been  developed  and  refined  as  a result.  The 
dependent  parameters  which  have  been  evaluated  can  be  summarized  under  four  general 
headings:  performance,  durability,  stability,  and  survivability.  The  five  design  variables  and  the 
four  dependent  variables  are  most  strongly  interrelated,  as  shown  in  the  matrix  of  Figure  90. 


Response  Variables 


Design  Variables 

Performance 

(Aerodynamic 

Efficiency) 

Durability 

(Bearing 

Life) 

Rotor 

Stability 

(Internal 

Friction) 

Survivability 
(Peak 
Transient 
Loads  and 
Deflections) 

a)  Support  System 

Stiffness 

X 

X 

X 

Damping 

X 

X 

X 

X 

Bumper  Clearance 

X 

X 

b)  Rotor 

Spline  Geometry 

X 

Spline  Lubrication 

X 

Figure  90.  Pertinent  Design  Factors  and  Their  Relevance  to  System  Response 


It  is  obvious  that  no  single  combination  of  design  variables  can  produce  the  best  response 
in  all  critical  areas  simultaneously.  For  example,  the  flexible  support  which  is  introduced  to 
decrease  dynamic  bearing  loads  may  dictate  increased  rotor/shr  oud  clearance  and  consequent 
performance  losses.  Thus,  the  “best”  design  involves  determining  the  most  beneficial  overall 
compromise.  It  is  also  obvious  that  the  qualities  that  are  most  desirable  are  unique  to  each  engine 
or  rig  because  of  its  specific  mission  or  function. 

It  is  evident  from  Figure  90  that  an  efficient  coordinated  design  system  is  required  to 
determine  the  “best  trade”  design.  This  system  was  developed  in  Task  X,  and  its  organization  is 
shown  in  Figure  91.  Inherent  in  this  system  is  the  recognition  that  the  term  "optimization” 
applies  to  bearing  life  and  aerodynamic  efficiency  effects.  It  is  also  recognized  that  these 
requirements  are  necessary  but  other  considerations  such  as  survivability  are  required  for  a truly 
optimized  design. 
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damped  support  parameters  to  achieve  maximum  bearing  life  and  to  minimize  rotor/shroud 
clearance  using  the  damper  analysis  developed  in  Task  I in  conjunction  with  steady-state 
dynamic  response  and  static  deflection  analyses,  and  (3)  verification  of  stability  using  the  Task 
IX  spline  analysis,  and  verification  of  survivability  using  a transient  response  analysis. 

To  establish  a starting  point  for  the  optimized  configuration,  the  following  procedure  is 
used.  The  rotor  dynamic  system  consisting  of  rotor,  bearings  and  bearing  supports  can  be  reduced 
to  an  approximate  single  mass  rotor  for  analysis  in  the  vicinity  of  the  first  shaft  bending  critical 
speed  (Reference  23).  This  is  accomplished  by  calculating  the  rigid  support  critical  speeds  of  the 
rotor  and  estimating  the  equivalent  rotor  stiffness  as  follows: 

K,  = M w*r  (142) 

where  M is  the  mass  inertia  of  the  rigid  rotor  and  wcr  is  the  first  rigid  support  critical  speed. 

Assuming  that  the  only  damping  in  the  system  comes  from  the  squeeze  film  damper,  the 
equations  of  motion  for  the  equivalent  single  mass  system  in  a flexible  damped  bearing  shown  in 
Figure  92  are 

MX,  + K,(X,  - X,)  = Mea>*e'“‘  (143) 

C.X,  + K,X[  + K,  (X,  - X.)  = 0 (144) 


23  Barrett,  L B..  E.  J.  Guntar,  Allair*.  ' Optimum  Bearing  and  Support  Damping  for  Unbalance  Rsaponaa  and  Stability 
of  Rotating  Mac  binary,"  ASMS  Paper  No.  77-GT-27. 
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Equations  (143)  and  (144)  are  combined  to  form  a single  equation  describing  the'  dynamic  motion 
of  the  mass: 


MX,  + C„  X,  + K*  X,  = MewV"*  (146) 

Cu  and  Km  are  the  effective  damping  and  stiffness  acting  on  the  rotor  mass.  The  expressions  for 
effective  damping  and  stiffness  are  as  follows: 


P _ K,»C, 

^ (K,+K,)«  + (wC,)* 

„ _ „ f K.fK.+K,)  + (tfC,) 

L (K,+K,):  + («C,)« 

Equations  (146),  (146),  and  (147)  can  be 
parameters: 


(148) 

(147) 


nondimensionalized  with  respect  to  the  following 


f = -^Ji- 
tter 


In  dimensionless  form,  Equations  (146),  (146),  (147)  become 

A 

X,  + 2 Wc,  X|  + w*  X,  ■=  eu  w*  e1-* 


(m  - 


A 

W* 


1> 

(1  + K)»  + (2f|,)' 

[ K (1  + K)  + (2fc,)«1 
L (1  + K)«  + (2fi,)*  J 


(148) 

(149) 

(150) 
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The  optimum  value  of  bearing  damping  ratio  maximizes  the  effective  damping  ratio,  £Je, 
and  the  optimum  bearing  damping  is  given  by 
t _ (1  ± K) 

£10  of  (151) 


and  the  maximum  effective  damping  ratio  is  given  by 


4f  (1  + K) 


Since  for  supercritical  speed  flexible  rotors  the  frequency  ratio  f is  unity,  the  expression  for 
optimum  damping  ratio  is 

1 + K 


and  for  maximum  effective  damping  ratio  is 


4 (1  + K) 


The  optimum  damping  ratio  is  plotted  in  Figure  93  for  determining  preliminary  optimum 
damping  coefficients. 


Figure  93.  Preliminary  Damper  Design  Curve 
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For  free  damped  vibrations  with  internal  viscous  friction  damping  such  as  spline  friction, 
Equation  (148)  becomes 


X,  + 2 u„  (£*m  + f,)  X,  + (w*  - 2i<i)crw  (,)  X = 0 
where  (,  is  the  internal  spline  friction  damping  ratio,  defined  as 


e.  = 


A. 


2Mwr 


(165) 


(156) 


C,  is  the  spline  coupling  damping  coefficient,  which  can  be  calculated  using  the  spline  friction 
model  developed  in  Task  VIII. 

The  solution  of  the  eigenvalues  of  Equation  (155)  reveals  that  the  rotor  system  will  only  be 
stable  if  the  operating  speed,  <*>,  is 


03  < 


Wcr 

(1  + K)' 


‘ 1 + 2K  > 

[i+  v1! 

L 2 J 

L 4f#  J 

(157) 


The  maximum  stable  speed  ratio,  ai/wrr,  at  which  the  rotor  can  operate  with  optimum  support 
damping  is  plotted  in  Figure  94  for  determining  preliminary  rotor  stability. 
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For  transient  blade  loss  excitation,  Equation  (148)  can  be  separated  into  an  x and  y 
direction  and  the  equations  can  be  solved  by  the  method  of  convolution  as  in  Reference  24. 
Equation  (148)  becomes 

X,  + 2 wcr  |j«  X»  + X,  = e.j  w*  coe  wt  (168) 

and 

Y,  + 2 wcr  Y,  + w*  Y,  - eu  «*  sin  wt  (169) 

The  maximum  transient  blade  loss  amplitude  is  predicted  by  calculating  the  maximum 
resultant  x,y  amplitude,  R,  by  the  use  of  the  convolution  integral: 

R = V X?  + Y/  (160  ) 

t 

v = -i-  f 0 at*  cos(wr)e  -(tMt-r)  ain<2(t— r)dr  (161) 

* *0  U 

I 

Y,  = cr  / euw»  sin(wr)e  -^(t'r)nnw<t-v)dr  (162) 

w o 

The  solutions  to  Equations  (160),  (161),  and  (162)  are  plotted  in  Figure  95  a8  the  ratio  of 
maximum  transient  amplitude  to  steady-state  amplitude  versus  speed  ratio.  The  curve  will  be 
used  for  determining  the  preliminary  rotor  transient  response  to  a sudden  blade  loss  unbalance. 

The  final  result  of  using  the  preliminary  design  curves  of  Figures  93  through  96  and  the 
design  optimization  method  of  Figure  91  will  be  the  definition  of  a power  turbine  rotor  dynamic 
system  which  (1)  has  the  best  possible  combination  of  bearing  life  and  aerodynamic  efficiency,  (2) 
is  stable  at  all  operating  speeds,  and  (3)  can  absorb  blade  loss  loads  without  damage  to  the 
engine. 
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Figure  95.  Preliminary  Transient  Response  Curve 


SECTION  VII 


CCNCLUSIONS  AND  RECOMMENDATIONS 


The  engine  rotor  dyne  nice  contract  contained  an  aggressive  analytical  and  complementary 
experimentel  program  to  develop  the  technology  required  for  designing  future  high  speed 
turboshaft  engines  with  optimum  rotor  dynamic  performancce.  The  design/prediction  system 
which  has  been  developed  can  be  used  to  design  future  gas  turbine  engine  rotors  to  have  improved 
efficiency,  durability,  and  survivability. 

The  squeeze  film  damper  characteristics  studied  were  analytically  predicted  and  ex- 
perimentally verified  to  be  strongly  dependent  on  damper  end-seal  leakage  and  inlet  feedback. 
The  damper  mathematical  model  developed  demonstrated  the  accuracy  end  versatility  necessary 
to  design  an  effective  squeeze  film  damper  with  any  combination  of  end  seal  and  inlet 
configurations  commonly  used  in  gas  turbine  engines.  However,  to  effectively  use  this  damper 
model,  additional  damper  tear  rig  data  is  required  which  will  quantify  the  leakage  characeristics 
and  parametrically  determine  the  best  type  of  end  seal  for  gas  turbine  engine  applications.  Also, 
a faster  solution  procedure  for  the  damper  model  should  be  pursued  in  order  to  improve 
computational  efficiency  for  transient  response  analyses  of  squeeze  film  damped  rotors. 

Mechanical  impedance  testing  of  three  types  of  flexible  bearing  supports  revealed  that 
satisfactory  accuracy  for  estimating  support  stiffness  can  be  accomplished  with  existing  beam 
theory  formulations.  If  a maximum  error  of  20%  cannot  be  tolerated,  then  a more  precise 
modeling  technique  such  as  a finite  element  method  should  be  used. 

The  synchronous  whirl  studies  verified  that  a properly  designed,  squeeze  film  damped, 
supercritical  rotor  can  have  acceptable  imbalance  sensitivity  and  can  sustain  transient  blade  loss 
imbalance  loads  without  damage  to  the  rotor  system.  The  existing  analytical  models  used  to 
predict  the  steady-state  and  transient  dynamic  response  of  synchronous  whirl  rig  gave  marginally 
satisfactory  results.  It  is  recommended  that  improved  computational  procedures  be  investigated 
in  order  to  reduce  the  numerical  errors  that  exist  in  commonly  used  engine  rotor  dynamic  models. 

The  nonsynchronouf  whirl  studies  verified  that  rotors  which  are  coupled  with  nonlubricated 
spline  couplings  are  prone  to  internal  friction  rotor  instability.  Lubricating  the  spline  coupling  or 
squeeze  film  damping  the  bearing  supports  has  a stabilizing  effect  on  the  rotor,  Furtheimore, 
squeeze  film  damping  can  reduce  the  steady-state  vibration  amplitude  at  the  critical  speed  by  as 
much  as  60%. 

While  it  is  felt  that  the  analytical  models  developed  under  this  program  along  with  existing 
rotor  dynamic  models  give  satisfactory  agreement  with  experimental  rig  data,  th“  analytical 
prediction  system  can  be  further  improved  with  additional  calibration  of  the  squeeze  film  damper 
with  empirical  damper  seal  leakage  data  from  parametric  rig  tests  of  commonly  used  seals.  Also, 
the  computational  precision  of  the  commonly  used  steady-state  and  transient  response  rotor 
dynamic  models  for  flexible  rotor  gas  turbine  engines  should  be  modified  to  use  improved 
numerical  methods  such  as  those  discussed  in  References  26  and  26. 


25  Hcrner,  G.  C , Pilkev  'The  Riccati  Transfer  Matrix  Method.”  ASME  Paper  No.  77-DET-32,  1977. 

26  Lund,  J.  W.,  “Modei  Responae  of  a Flexible  Rotor  ir.  Fluid-Film  Bearings,"  ASME.  Journal  of  Engineering  for 
Industry,  page  625,  f.i  ..y  1974. 
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APPENDIX  A 

DAMPER  ANALYSIS  BY  NUMERICAL  DIFFERENTIATION  AND  INTEGRATION 

1.0  GENERAL  INFORMATION 

This  document  contains  the  operating  instructions  for  the  Damper  Analysis  by  Numerical 
Differentiation  and  Integration  (DANDI)  Computer  Program  developed  under  Contract 
DAAJ02-76-C-0011,  Engine  Rotor  Dynamics.  Descriptions  of  input  and  output  program  data  are 
provided  along  with  flow  diagrams. 

Program:  Damper  Analysis  bv  Numerical  Differentiation  and  Integration 
(DANDI) 

Language:  FORTRAN  IV/Level  221 
Computer  Configuration:  IBM  370/168 
Operating  System:  OS  370,  Release  3.7VS 
Estimated  Running  Time:  0.2  sec  CPU  time 

Program  Author:  Ronald  A.  Marmol,  Pratt  & Whitney  Aircraft  Group, 

P.  O.  Box  2691,  West  Palm  Beach,  Florida  33402 

2.0  INPUT  DESCRIPTION 

Figure  A-l  is  a sample  data  input  form  showing  the  data  input  format.  Table  A-l  lists  the 
input  data  definition. 

3.0  OUTPUT  DESCRIPTION 

Table  A-2  lists  the  output  data  definition.  Figure  A-2  is  a flow  chart  for  the  computer 
program  and  Figure  A-3  is  a logic  flow  diagram  for  choosing  boundary  condition  subroutines. 


TABLE  A-l 

INPUT  DATA  DEFINITION 

Card 

No.  Cols.  Type  Description 


1-80 

A.N. 

Any  descriptive  title  information. 

1-10 

F.P. 

RB  = 

Radius  of  bearing  (in.). 

11-20 

F.P. 

C - 

Damper  bearing  clearance  (in.). 

21-30 

F.P. 

EL  = 

Damper  bearing  length  (in.). 

31-40 

F.P. 

ELI  = 

Axial  dimension  from  the  inlet  plane  to  the  closest 
end  (in.).  If  the  damper  has  symmetry,  it  must  have 
the  inlet  centered  between  identical  end  conditions. 

41-44 

Int. 

NI  - 

Type  of  inlet  (right  adjusted  with  no  decimal  point). 

0 or  left  blank  for  no  inlet. 

1 for  single  inlet  hole  with  backflow. 

2 for  two  inlet  holes  180  deg  apart. 

• 

£ 

• 

• 

100  for  circumferential  groove. 

46 

Int. 

NSl  = 

s 

Type  of  end  seal  at  Z=0.0  (right  adjusted  with  no 
decimal  point). 

1 for  piston  ring  seal. 

2 for  radial  O-ring  seal. 

3 for  side  O-ring  seal. 

4 for  no  seal. 

48 

Int. 

NS2  - 

Type  of  end  seal  at  Z = L (right  adjusted  with  no 
decimal  point).  Same  options  as  NSl  above. 

60 

Int. 

NSYMM  = 

s 

as 

Symmetry  option  (right  adjusted  with  no  decimal 
point). 

0 or  left  blank  for  no  symmetry. 

1 for  inlet  centered  between  two  identical  end  seals 
(NI  * 0). 

52 

Int. 

ITERM  * 

sx 

Simplified  equation  option  (right  adjusted  with  no 
decimal  point). 

0 or  left  blank  for  complete  equations  to  be  used. 

1 for  simplified  equations  to  be  used  for  FI,  F3  and  F4. 

1-16 

F.P. 

EMU  = 

Lubricant  viscosity  (lb -sec/in. **2). 

16-30 

F.P. 

RPM  - 

Shaft  spaed  (rpm). 

31-46 

F.P. 

PSUPLY  - 

Supply  pressure  (pei).  Mint  be  non-zero. 

46-60 

F.P. 

PEXIT  - 

Exit  pressure  (pei). 

61-76 

F.P. 

PCAV  - 

Cavitation  pressure  (pei). 
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TABLE  A-l 

INPUT  DATA  DEFINITION  (Continued) 


Card 

No.  Cols.  Type 

4.  1-4  Int. 

5-8  Int. 


9-12  Int. 


Description 

M = Number  of  circumferential  grid  points.  Right  ad- 
justed with  no  decimal  point.  Even  Number  — 
Maximum  of  60. 

N = Number  of  axial  grid  points.  Right  adjusted  with  no 
decimal  point. 

Odd  number. 

If  NSYMM  = 0,  Use  5 £ N £ 15. 

If  NSYMM  = 1,  Use  N = 9,  13,  17,  21,  25,  or  29. 

I0P  = Solution  option  (right  adjusted  with  no  decimal 
point). 

= 1 for  transient  damper  forces  (FX  and  FY). 

= 2 for  linearized  squeeze  film  stiffness  (K),  squeeze 
film  damping  (C),  and  damper  offset  (EP). 

= 3 for  Tondl  coefficients  (Fl,  F2,  F3,  and  F4). 

= 4 for  Lund  damoing  coefficients  (CXX,  CXY,  CYX, 
and  CYY). 

- 5 for  transient  damper  forces  (FX  and  FY)  and  grid 
point,  pressures. 


If  IOP  = 4,  input  the  following  card  No.  5. 


6. 

1-15 

F.P. 

*X 

■= 

Dimensionless  journal  position  (X/C). 

16-30 

F.P. 

«Y 

= 

Dimensionless  journal  position  (Y/C). 

31-46 

F.P. 

AcX 

ae 

Dimensionless  journal  velocity  (VX/WC) . 

46-60 

F.P. 

A«Y 

- 

Dimensionless  journal  velocity  (VY/WC). 

6. 

1-15 

F.P. 

CP 

= 

"Ix>ss  coefficient”  for  flow  through  the  radial  gap 
around  the  piston  ring  seal.  0.0  £ CP  ^ 1.0 

16-30 

F.P. 

HS 

St 

Radial  gap  around  the  piston  ring  seal  (in.). 

31-45 

F.P. 

WP 

Width  (axial  dimension)  of  the  piston  ring  (in.). 

7. 

1-15 

F.P. 

COR 

a; 

Flow  coefficient  for  O-ring  seal  (in  *J3/8ec)  / (psi-in). 

16-30 

F.P. 

HOR 

= 

Radial  thickness  of  O-ring  as  installed  in  groove  (in.). 

31-45 

F.P. 

CKOR 

- 

O-ring  stiffness  — axial  defl.  / psi.  (in./psi). 

8. 

1-15 

F.P. 

Cl 

Flow  coefficient  for  one  inlet  hole  (in.**3/sec)  / psi. 

TABLE  A-2 

OUTPUT  DATA  DEFINITION 


Printout Description 

(Title)  Any  descriptive  title  information. 

Bearing  Radius  Radius  of  bearing  (in.). 

Inches 

Damper  Bearing  Damper  bearing  clearance  (in.). 

Clearance 

Inches 

Axial  Dimension  From  Axial  dimension  from  the  inlet  plane  to  the  closest  end  (in.).  If  the 

Inlet  Plane  to  Closest  damper  has  symmetry,  it  must  have  the  inlet  centered  between 

End  Inches  identical  end  conditions. 


Inlet  Type  of  inlet 

Type  0 or  left  blank  for  no  inlet. 

1 for  single  inlet  hole  with  backflow. 

2 for  two  inlet  holes  180  deg  apart. 

• 


Z=0 

End 

Seal 

Type 


Z = L 
End 
Seal 

Type 

Symmetry 

Option 


Equation 

Option 

Lubricant 
Viscosity 
ib -sec/in.  **2 

Shaft  Speed 
rpm 

Supply 

Pressure 

psi 


100  for  circumferential  groove. 

Type  of  end  seal  at  Z=0.0 

1 for  piston  ring  seal. 

2 for  radial  O-ring  seal. 

3 for  side  O-ring  seal. 

4 for  no  seal. 

Type  of  end  seal  at  Z = L.  Same  options  as  NS1  above. 


Symmetry  option. 

0 or  left  blank  for  no  symmetry. 

1 for  inlet  centered  between  two  identical  end  seals  (HI  * 0). 

Simplified  equation  option. 

0 or  left  blank  for  complete  equations  to  be  used. 

1 for  simplified  equations  to  be  used  for  FI,  F3  and  F4. 

Lubricant  viscosity  (lb-sec/in.**2). 


Shaft  speed  (rpm). 

Supply  pressure  ipsi).  Must  be  ion-*ero. 


134 


TABLE  A-2 

OUTPUT  DATA  DEFINITION  (Continued) 


Printout 

Description 

Exit 

Pressure 

Exit  pressure  (psi). 

psi 

Cavitation 

Pressure 

Cavitation  pressure  (psi). 

psi 

Number  of 
Circumferential 
Grid  Points 

Number  of  circumferential  grid  points.  Even  number  — maximum  of 
60. 

Number  of 

Axial 

Grid  Points 

Number  of  axial  grid  points.  Odd  number. 

If  NSYMM  = 0,  Use  5 < N < 15. 

If  NSYMM  = 1,  Use  N = 9,  13,  17,  21  25,  or  29. 

Solution 

Option 

1 for  transient  damper  forces  (FX  and  FY) . 

2 for  linearized  squeeze  film  stiffness  (K),  squeeze  film  damping  (C), 
and  damper  offset  (EP). 

3 for  Totidl  coefficients  (FI,  F2,  F3,  and  F4). 

4 for  Lund  damping  coefficents  (CXX,  CXY,  CYX,  and  CYY). 

5 for  transient  damper  forces  (FX  and  FY)  and  grid  point  pressures. 

If  Solution  Option  = 1 or  5 

Journal  Position 

X 

X position  of  the  damper  journal  center  (in  ). 

Y 

Y position  of  the  damper  journal  center  (in.). 

Journal  Velocity  (inches/sec) 

VX 

Velocity  of  the  damper  journal  center  in  the  X direction  (in./sec). 

VY 

Velocity  of  the  damper  journal  center  in  the  Y direction  (in. /sec). 

If  Solution  Option  = 2 

Damper 

Static  Load 
lb 

Damper  static  load  (lb). 

Whirl 

Eccentricity 

Irches 

Whirl  eccentricity  (in.). 

Mechanical 
Spring 
Stiffness 
lb /in. 

Mechanics'  spring  stiffness  (lb /in.). 
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TABLE  A-2 

OUTPUT  DATA  DEFINITION  (Continued) 

Printout Description 

If  Solution  Option  = 3 

Radial  Radial  velocity  (in./sec). 

Velocity 

Inches/sec 

Damper  Damper  eccentricity  (in.). 

Eccentricity 

Inches 


If  Solution  Option  = 4 


Dimensionless  Journal  Position 

Y/C  Dimensionless  journal  position  (X/C). 

Y/C  Dimensionless  journal  position  (Y/C). 

Dimensionless  Journal  Velocity 

VX/WC  Dimensionless  journal  velocity  (VX/WC). 

VY/WC  Dimensionless  journal  velocity  (VY'/WC). 

For  Solution  Option  = 1,  2,  3,  4,  or  5 


Piston  Ring 
Seal  Loss 
Coefficient 

Piston  Ring 
Seal 

Radial  Gap 
Inches 

Piston  Ring 
Width 
Inches 

O-Ring 
Seal  Flow 
Coefficient 
(in.**3/sec)/(psi-in.) 


‘Loss  coefficient’  for  flow  through  the  radial  gap  around  the  piston  ring 
seal.  0.0  <;  CP  £ 1.0 


Radial  gap  around  the  piston  ring  sea)  (in.). 


Width  (axial  dimension)  of  the  piston  ring  (in.). 


Flow  coefficient  for  O-ring  seal  .'in.**3/sec)/(psi-in.). 
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Printout 

0-Ring 

Radial 

Thickness 

Inches 

0-Ring 

Stiffness 

in./psi 

Inlet 

Hole  Flow 
Coefficient 
(in,**  3/sec  i/psi 

EPS 

PHI 

EDOT 

PHIDOT 

DEPSDT 

DPHIDT 

PS 

PE 

OMEGA 


FX 

FY 


Squeeze 
Film  Stiffness 
lb  /in. 

Squeeze 
Film  Damp 
lb-sec/in. 

Damper 

Offset 

Inches 


TABLE  A-2 

OUTPUT  DATA  DEFINITION  (Continued) 

Description  

Radi"!  thickness  of  O-ring  as  installed  in  groove  (in.). 


O-ring  stiffness  = axial  defl./psi  (in./psi). 


Flow  coefficient  for  one  inlet  hole  (in.**3/sec)/psi. 


(whirl  eccentricity)/(damper  bearing  clearance) 

Damper  angular  position  in  clearance  circle  (radians). 

Radial  velocity  of  damper. 

Whirl  velocity  of  damper. 

(Radial  velocity  of  damper)/(bearing  clearance  X shaft  angular 
velocity) 

(Whirl  velocity  of  damper)/(shaft  angular  velocity) 

(Supply  pressure)  *C**2/(6.0  * RMU  * OMEGA  * RB  **2) 

(Exit  pressure) *C * *2/(6.0  * EMU  * OMEGA  * RB  **2) 

(Shaft  rpm)  * ir/30. 

If  Solution  Option  = 1 

Resulting  damper  force  in  the  horizontal  direction. 

Resulting  damper  force  in  the  vertical  direction. 

If  Solution  Option  - 2 

Damper  Btiffness  coefficient  for  circular  whirl. 

Damper  damping  coefficient  for  circular  whirl. 

Damper  offset  from  center  of  clearance  circle  due  to  static  loading. 
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TABLE  A -2 

OUTPUT  DATA  DEFINITION  (Continued) 

Printout Description 

If  Solution  Option  = 3 

FI  Tondl  coefficients  as  defined  by  Reference  2. 

F2 

F3 

F4 

If  Solution  Option  = 4 


CXX  Linear  damping  coefficients  as  defined  by  Lund 

CXY 

CYX 

CYY 

If  Solution  Option  = 5 

Finite  difference  grid  point  pressures  are  printed 
FX  Resulting  damper  force  in  the  horizontal  direction. 

FY  Resulting  damper  force  in  the  vertical  direction. 


4 
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(J  M ^b’  L,  L-),  Pg,  Pg,  C| 

M.  w,  Xj,  Yj,  Xj,  Vj,  Cp,  hs,  Cqr.  Kqr.  hoR 
Inlet  No.  "N1",  Seal  No.  NS1,  NS2 
NSYMM  (Axial  Symmetry) 


Calculate  *■  Groups,  Position, 
Motion,  Digitize  i (J),  H(J) 

_d  d£  C C R) 

*•  dT'  0,dT  Rj*  Rb‘  Rb 


Calculate  All  Matrix  Elements 
(Except  End  Boundary  Rows) 
from  Reynold's  Equation, 
Ajk,  and  Right  Side  Vector 


Call  Proper  Subroutines 


Calculate  Matrix  Elements 
for  Boundaries  at  Each  End 
(Seals),  and  at  Inlet 


Return 


Print  Out  Matrix  Elemen*s 
Ajk  and  Solution  Vector  P JK 


- 1 * 


APPENDIX  B 

SPLINE  COUPLING  ANALYSIS  PROGRAM 

1.0  GENERAL  INFORMATION 

Program:  Spline  Coupling  Analysis  Program 

Language:  FORTRAN  IV/Level  221 
Computer  Configuration:  IBM  360/370 
Operating  System:  SCOPE  3.4 
Estimated  Running  Time:  0.5  sec  CPU  time 

Program  Author:  Judith  Hutchison,  Joseph  A.  Tecxa,  Mechanical  Technology  Inc. 
968  Albany  Shaker  Rd.,  Latham,  N.Y.  12110 

2.0  INPUT  DESCRIPTION 

The  computer  software  developed  for  this  project  consists  of  a main  program  and  five 
subroutines:  SPLININ,  SPLINE,  SPLOUT,  OUTLINE,  and  QB.  The  subroutines  are  designed 
to  be  inserted  into  a rotor  dynamics  program  which  is  simulated  by  the  short  main  program.  The 
total  input  (exclusive  of  control  cards)  consists  of  seven  cards  — three  read  by  the  main  program 
and  the  remaining  four  by  subroutine  SPLININ.  These  cards  are  described  in  the  subsections 
which  follow. 

2.1  Main  Program  Input 

Program  MAIN,  which  calls  SPLININ  and  SPLINE,  requires  three  inputs,  supplied  by  3 
cards.  The  quantities  in  parentheses  indicate  the  card  formats. 

Card  1 — Title  Card 

TITLE  (I),  1=1,8  (8  A10) 

Card  2 — Type  of  Fit 

IFTT  = 1 for  major  diameter  fit  (1615) 

= 2 for  side  fit 

Card  3 — Operating  Speed 

SPEED  (rpm)  (8E10.3) 

2.2  Subroutine  SPLININ  Input 

Card  l — Integer  Controls  (1615) 

NT  — the  number  of  tooth  pairs  in  the  coupling. 

NK1  — the  number  of  sections  into  which  the  external  tooth  is  divided  for 
analysis  (must  be  £25;  10  recommended). 
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NK2  — the  number  of  sections  into  which  the  internal  tooth  is  divided  for 
analysis  (must  be  £25;  10  recommended). 

Many  of  the  dimensions  referred  to  on  the  following  cards  are  described  in  the  standard 
INVOLUTE  SPLINES  AND  INSPECTION,  ANSI  B92.1,  1970,  English  Version,  Published  by 
Society  of  Automotive  Engineers,  Inc.,  Two  Pennsylvania  Plaza,  N.Y.,  N.Y.  10001. 

Card  2 (8E10.3) 

PSI  — pressure  angle  at  pitch  radius  (degrees). 

XL  — face  width  (inches).  See  figure  B-l  for  detail  view. 


t 

I 


Figure  B-l.  Views  of  Spline  Coupling 
EMOD  — Young’s  modulus  (lb/sq  in.). 

XNU  — Poisson’s  ratio. 

WC  — width  of  contact  patch  (inches).  This  is  the  radial  distance  along 
which  contact  occurs  between  male  and  female  teeth.  For 
involute  spline  teeth  it  is  the  difference  in  form  radii  of  the 
external  and  internal  spline  teeth. 

EPS  — assumed  angular  whirl  amplitude  of  coupling  (degrees).  This  is 
the  relative  angle  between  internal  and  external  coupling  halves 
of  the  coupling  during  unstable  whirl.  This  value  must  be 
assumed  and  it  is  expected  to  be  slightly  larger  than  the  angle 
defined  by  half  the  face  width  and  the  radial  clearance. 

EPS  £ tan-1  [ 2 y.gkgAR.] 
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OMEGAl  — first  natural  frequency  (rpm).  Unstable  whirl  is  expected  to 
occur  at  the  first  natural  frequency.  OMEGAl  must  be  esti- 
mated and  input  to  allow  correct  calculation  of  equivalent 
viscous  damping  coefficients. 

DELTAO  — assumed  radial  whirl  amplitude  (inches).  This  is  t,he  relative 
radial  amplitude  between  internal  and  external  spline  members 
expected  during  unstable  whirl.  This  assumed  value  is  expected 
to  be  incrementally  larger  than  the  radial  clearance. 


Card  3 (8E10.3) 


XMU  — coefficient  of  friction.  This  value  may  vary  from  approximately 
0.06  to  0.26,  depending  on  the  type  of  lubrication  and  the  relative 
sliding  velocity. 


TORQ  - 
CLEAR  - 
G - 
S — 

T - 


transmitted  torque  (in.- lb) 
radial  clearance  (in.- lb) 
shear  modulus  (lb/sq  in.) 
space  width  (in.) 

tooth  thickness  (in.) 


These  are  described  in 

} ANSI  B92.1 1970,  the  stan- 
dard referred  to  earlier. 


Card  4 (8E10.3) 

All  values  on  this  card,  except  for  the  fillet  radii,  may  be  taken  from  ANSI  B92.1 
1970  and  are  shown  schematically  in  figure  B-2. 


D 

DRE 

DFE 

RFILl 

DRI 

DFI 

RFIL2 


— pitch  diameter  (in.) 

— minor  diameter,  external  (in.) 

— form  diameter,  external  (in.) 

— fillet  radius,  external  (in.) 

— major  diameter,  internal  (in.) 

— form  diameter,  internal  (in.) 

— fillet  radius,  internal  (in.) 


Figure  B-3  shows  a coding  sheet  which  describes  the  location  of  each  input  quantity  in  the 
card  deck. 


9.0  OUTPUT  DESCRIPTION 

Page  1 — Complete  list  of  input  parameters 


Page  2 — (Numbers  in  parentheses  refer  to  equations  in  spline  coupling 
stability  report,  which  is  the  reference  for  this  program.) 


a.  Radial  stiffness,  angular  stiffness,  radial  damping, 
angular  damping 

b.  Contracted  stiffness  matrix 


c.  Contracted  damping  matrix 

d.  Expanded  stiffness  matrix  (2-28) 


} 


(2  26) 


e.  Expanded  damping  matrix  (2-29) 

Figures  B-4  and  B-5  contain  flow  charts  of  the  main  program  and  the  spline  coupling 
subroutine. 


A.  Major  Diameter  Fit 


Figure  B-2.  Input  Values  of  Card  4 


Figure  B-3.  Spline  Program  Coding  Sheet 


Spline 


Outline 

Calculate  Tooth 
Half  Widths  as 
Functions  of  Root 
and  Fillet  Radii 


For  Male  Tocth 
Calculate: 

1)  A'ea  of  Slice 

2)  Second  Moment  of 
Area 

3)  Moment  Arm 

TL 

For  Female  Tooth 
Calculate: 

1)  Area  of  Slice 

2)  Second  Moment  of 
Area 

3)  Moment  Arm 

Mojor  Diameter  \ Yes 
Fit?  / H 


Calculate  Base 
Rotation  Compliance 


Compliance  Due 
to  Tooth  Bending 


Calculate 
Tooth  Stiffness 


Calculate  Compliances 
Due  to: 

1)  Shear 

2)  Base  Rotation 

3)  Contact 
Deformation 


Calculate 
Lateral  Stiffness 


Calculate 
Tooth  Stiffness 


Calculate 
Lateral  Stiffness 


£]£»] 


4.0  EXAMPLE  CASES 


Example  cases  have  been  run  for  one  of  each  coupling  type  — major  diameter  fit  and  side 
fit.  Table  B-l  lists  the  input  quantities  for  each  of  the  two  cases;  the  results  for  the  example  are 
given  in  table  B-2. 


Table  B-l.  Input  Data  for  Two  Example  Cases 


Data  For 
Major  Diameter 
Fit  Coupling 

Data  For  Side 
Fit  Coupling 

Operating  Speed  (rpm) 

4800 

4900 

No.  of  Tooth  Pair* 

36 

36 

No.  at  Segments,  Ext. 

10 

10 

No.  at  Segments,  Int. 

10 

10 

Pressure  Angle  (deg) 

30 

30 

Face  Width  (in.) 

1.76 

1.76 

Young’s  Modulus  (lb/aq  in.) 

30X10* 

30X10* 

Poisson's  Ratio 

0.3 

0.3 

Contact  Patch  Width  (in.) 

0.1 

0.1 

Angie  of  Tilt  (deg) 

0.1026 

0.3448 

Pint  Critical  (rpra) 

2400 

2400 

Radial  Amplitude  (in.) 

1.662X10"' 

6.262X10"* 

Coefficient  of  Friction 

0.16 

0.15 

Torque  (in. -lb) 

2000 

2000 

Radial  Clearance  (in.) 

1.66X10"* 

5.26X10"* 

Shear  Modulus  (lb/sq  in.) 

11.68X10* 

11.68X10-* 

Space  Width  (in.) 

0.1344 

0.1344 

Tooth  Thickness  (ext)  (in.) 

0.1268 

0.1239 

Pitch  Dia.  (in.) 

3.0 

3.0 

Minor  Dia.  (ext)  (in.) 

2.871 

2.83 

Form  Dia.  (ext)  (in.) 

2.910 

2.91 

Fillet  Radius  (ext)  (in.) 

0.040 

0.040 

Major  Dia.  (int)  (in.) 

3.084 

3.170 

Form  Dia.  (int)  (in.) 

3.009 

3.090 

Fillet  Radius  (int)  (in.) 

0.040 

0.040 

Table  B-2.  Results  for  Two  Example  Cases 

Major  Diameter 
fit  Output 

Side  Fit  Output 

Radial  Stiffness  (IbAn.) 

1.2978X10* 

1.3161X10* 

Angular  Stifhsaa  (in.  lb/radian) 

3.2970X10* 

S.SHSXIO1 

Radial  Damping  (lb- sac/in.) 

417.82 

133.47 

Angular  Damping  (in.lb-eec/radian) 

168.68 

47.143 
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